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Abstract. Let (SDq) be the Stokes operator defined in a bounded domain Q of M 3 with 
Dirichlet boundary conditions. We prove that, generically with respect to the domain Q with 
C 5 boundary, the spectrum of (SDq) satisfies a non resonant property introduced by C. Foias 
and J. C. Saut in [TB] to linearize the Navier-Stokes system in a bounded domain Q of R 3 
with Dirichlet boundary conditions. For that purpose, we first prove that, generically with 
respect to the domain Q with C 5 boundary, all the eigenvalues of (SDq) are simple. That 
Q_| ! answers positively a question raised by J. H. Ortega and E. Zuazua in [251 Section 6]. The 
<^ \ proofs of these results follow a standard strategy based on a contradiction argument requiring 
r-| ; shape differentiation. One needs to shape differentiate at least twice the initial problem in the 
direction of carefully chosen domain variations. The main step of the contradiction argument 
amounts to study the evaluation of Dirichlet-to-Neumann operators associated to these domain 
variations. 
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1 Introduction and main results 

In this paper, we consider the eigenvalue problem for the Stokes system with Dirichlet bound- 
ary conditions defined in a bounded open subset Q of M 3 with C e boundary, £ > 4, 



(SD S 



-A0 + Vp 
div</> 



p(x)dx 



\cf) in Q, 

in Q, 

on <9f2, 

0. 



Here we use <fi G M. 3 and p G K to denote respectively the velocity field and the pressure. It is 
well-known that (SD^) admits an increasing sequence of positive eigenvalues (A n ) n >i tending 
to infinity as n goes to infinity. 

The purpose of this paper is to prove genericity results on the spectrum of (SD^) with 
respect to domains of R 3 . We start by clarifying the notion of genericity considered below. 
Recall that the set of bounded domains of M 3 with C boundary denoted by D 3 can be endowed 
with the following topology: the base of open neighborhoods is (essentially) given by the sets 
V(Q, e) defined, for any domain Q G Of and e > small enough, as the images of Q by Id + u, 
with u G W e+1 '°°(n,R 3 ) and HuH^+i,. < e (cf. [20] and [29]). Then e is chosen so that 
Id + u : Q — > (Id + u)(Jl) is a diffeomorphism. As shown by A. M. Michelleti in [23] (see also 
[2"Ul Appendix 2]), a neighborhood V(Q, e) of Q G D| as defined previously is metrizable using 
a Courant-type distance, denoted by o^+i, and each (V(Q,e),di+i) is complete and separable. 
For any domain Q G D 3 , we use I0 3 (fi) to denote the Banach manifold obtained as the set of 
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images (Id + u)(Q) by u G W e+1,oc (Q, IR 3 ), which are diffeomorphic to Q. A property (P) will 
thus be referred to as "being generic with respect to Q G D 3 " if, for every Q G Of, the set of 
of domains of Df (f2) where (P) holds true contains a countable intersection of open and dense 
subsets of Df(fi). 

The main contribution of this paper consists in proving a conjecture formulated by C. 
Foias and J. C. Saut in [IB] on the generic non-resonant character of the spectrum of the 
Stokes operator with Dirichlet boundary conditions. Note that the non-resonance condition 
plays a crucial role in linearizing the Navier-Stokes systems to obtain a normal form, and 
then computing a useful asymptotic expansion for its solutions (cf. [16], [12] and [13]). This 
situation is analogous to the classical Poincare's normal form theorem for ordinary differential 
equations (cf. [5J Chapter 5]) although the proof is more involved. See also [H] for a recent 
development in the spirit of [IB] . As noticed in [16], non-resonance does not occur for periodic 
boundary conditions. However, the authors conjectured that non-resonance should be generic 
for Dirichlet boundary conditions. In this paper, we confirm that conjecture in Theorem 11.11 
Let us first recall the following definition (cf. [161 Definition 1]). 

Definition 1.1. We call resonance in the spectrum of (SDq) a relation of the type 



If no resonance occurs in the spectrum of (SDq), then (SDq) will be called non-resonant. 

Theorem 1.1. Generically with respect to Q G Df, the spectrum of the operator (SDq) is 
non-resonant. 

In order to prove Theorem ll.l| we need first to establish another genericity result on the 
spectrum of the Stokes operator known as generic simplicity. 

Theorem 1.2. Generically with respect to Q G D|, all the eigenvalues of (SDq) are simple. 

Remark 1.1. In [25], several properties for the Stokes system with Dirichlet boundary con- 
ditions (in particular the simplicity of the spectrum) were proved to be generic for domains 
in M 2 . Moreover, in the same paper, the three dimensional case was considered in Section 6, 
pointing out why techniques developed in [25] could only handle the two dimensional case. In 
this regard, Theorem 11.21 answers positively the open question of Section 6 in [25]. 

Theorem 11.21 is of course a particular case of Theorem ll.l| but will allow us to work only 
with simple eigenvalues in the main step of Theorem II . If s proof. This reduction is essential 
in our arguments. We now describe the strategy of the proofs. As it is standard since [I], the 
reasoning goes by contradiction and is based on shape differentiation. 

We start with a description of the proof of Theorem 11.21 Fix a domain Qq G D 3 . For every 
integer k, we define as the (open) subset of D 3 (f2 ) whose elements Vt verify that the k first 
eigenvalues of (SDq) are simple. Clearly, by Baire's lemma, proving Theorem 11.21 amounts to 
show that Ak+i is dense in Ak for every k > 0. We argue by contradiction and assume that 
there exists an integer k, a domain Q with C boundary in Ak and e > 0, such that, for every 
u G W e+1,oc (Q, IR 3 ) with HiiHtyf+1,00 < e, the domain (ld + u)Q, or simply Q + u, belongs to Ak 
but not to A k+ i- Let m > 2 be the multiplicity of A, the value of the (k + l)-th eigenvalue of 
(SDq) and 0j, i — 1, ■ ■ • , m, orthonormal eigenfunctions associated to A. Finally, let n be the 



k 




where rrij G N, 



1 < j < k. 



(1) 
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outward unit normal vector field on dQ. By computing the shape derivative of the (n + l)-th 
eigenvalue of (SDq), J. H. Ortega and E. Zuazua obtained in [25] that, at every x G dQ, one 
has, for i,j — 1, • • ■ , rn, and i ^ j, 

^ti. n = o ^ti.^tl = o II — f = II — f (2) 
dn dn dn dn dn 

dcj) 

If m > 2, then there necessarily exists 1 < i < m so that -—— = on dQ and one reaches a 

dn 

contradiction using a unique continuation result due to Osses (cf. [27]). However, in order to 
obtain generic simplicity (rn — 1), it was not clear how to pursue the reasoning by contradic- 
tion, i.e., showing that relations in ([2]) do not hold true generically with respect to the domains 
of R 3 if m = 2. Note that, for questions involving scalar PDEs, if one wants to prove generic 
simplicity of the spectrum of a self-adjoint operator with Dirichlet boundary conditions, then 
it is standard to follow the lines of the above mentioned contradiction argument and to reach 
Eq. ([2]). The second equation there is now a product of real numbers and a contradiction 
follows readily by unique continuation, cf. [I] and [20]. Therefore, the difficulty for showing 
the generic simplicity of the spectrum of (SDq) stems, at this stage of the argument, from the 
vectorial character of 0,, i.e., the fact that we are dealing with a system of PDEs. 

In this paper, we push further the contradiction argument by computing the shape deriva- 
tive of the (k + l)-th eigenvalue of (S)q +u at every u G W £+1 '°°(Q, R 3 ) with || 
small enough. The relations obtained in Eq. fl2]) for Q are now valid for every domain Q + u 
with u small enough. At this stage, we are not able to derive a contradiction. So we again 
take the shape derivative of the above relations on dQ and end up with expressions of the type 

M\u)(x) = -(u-n)(x)^^(x), xedQ, (3) 



for w e+i,co < e and where 



d(j)i d4>j M<90 i||2 ,.<% 2 



^ dn ' dn dn ' ^ dn ^ 1 dn " 

Taking into account the expression of M, its shape derivative M'[u) can also be expressed 
in terms of Neumann data of the shape derivatives of the eigenfunctions whose values on 
dQ have the regularity of u • n. By standard elliptic theory, if u ■ n belongs to the Sobolev 
space H s (dQ), M'(u) a priori belongs to H s ~ 1 (dQ). Then, the key observation is that a 
gap of regularity exists between the two sides of Eq. ([3]) since the right-hand side trivially 
belongs to H s (dQ), for s < 1 and £ > 4, the latter assumption needed to assert that 9R g^ 
is continuous on dQ. The whole point now comes down to use that gap of regularity in 
order to reach a contradiction. In this paper, we reformulate the issue at hand as follows: 
how to extract pointwise information (i.e., for x G dQ) reflecting the aforementioned gap of 
regularity and thus allow us to pursue the reasoning by contradiction. This rather elementary 
line of attack, first considered in [S] and also applied in [B], consists in choosing appropriate 
variations u "localized" at an arbitrary point x G dQ. We note that problems treated in 
[9] and [6J concerned planar domains and, therefore, equations of the type ([3]) were valid on 
closed C 3 curves of R 2 . In that case, the localization procedure is easier to handle. Indeed, 
the strategy adopted in [9] and [6J consisted in extending M'{u) for variations u defined on 
dQ as continuous functions except at some point x G dQ. More particularly, u = u x can be 
taken as a Heaviside like function admitting a single jump of discontinuity at x. In order to 
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exploit the gap of regularity, the singular part of M'(u x )(-) at x (in the distributional sense) 
had to be computed, to eventually obtain the following expression, 



where a denotes the arclength (with a = corresponding to x) and R(-) belongs to H 1 ^ 2 ^ e (dVt) 
for every e > 0. Plugging back the above expression into Eq. (jSJ), one deduces that M (-) = 
on dfl. In [9], the previous relation provided additional information and allowed to conclude 
the contradiction argument. However, in [Hj, it turns out that M (-) is proportional to M(0)(-) 
and hence is trivially equal to zero. To determine the first non trivial term in the "singular" 
expansion of M'{u x ) + (u x ■ n) 9A q^ at x, in the distributional sense, a detailed study of 
Dirichlet-to-Neumann operators associated to several Helmholtz equations was required. 

In the present paper, the "localization" procedure, i.e., the choice of appropriate variations 
u for any arbitrary point x G dQ, must be performed for functions defined on a surface dQ 
and not anymore on a curve, as in [Hj and [5J. For that purpose, after fixing an arbitrary 
point x G dQ, we will choose sequences of smooth functions u £ ^ Xe approximating the Dirac 
distribution at x £ as e tends to zero, the point x £ G dQ being any point at distance e of x. The 

gap of regularity between the two sides of Eq. Q will be now quantified in terms of powers of - 

e 

and not anymore measured in the distributional sense. We are therefore led to an asymptotic 
analysis as - tends to infinity and, more precisely, the right-hand side of Eq. ([3]) is a O(-) 

w 

meanwhile we will establish that the left-hand side of of Eq. ([3]) is equal to — ^ + 0( — ), where 

e 6 e 2 

w £ is bounded independently of e. Letting e tend to zero, one deduces that lim^o w e = 
and finally one concludes the contradiction argument and Theorem 11.21 is established. Note 
that the exact characterization of w £ requires, as in [5J, a detailed study of certain Dirichlet- 
to-Neumann operators, but here, associated to the Stokes system. That study heavily uses 
many technical results borrowed from [201 Chap. 7], not only for handling certain weakly 
singular operators but also for the material which is necessary to evaluate integrals defined 
on the surface dVt. It is noteworthy that, to perform the evaluation of the surface integrals, 
we choose charts based at x £ G dVt near the fixed point x G dQ, but not exactly at x. This 
trick turns out to be crucial for handling the singularities in computations involving boundary 
layer potentials. Of an equal importance, it also provides two degrees of freedom, namely the 
distance and the angle (in local coordinates) between x £ and x, and functions of these two 
variables being equal to zero give additional information to yield a contradiction. 

Let us now briefly mention how the argument for Theorem 11.11 goes. Since the resonance 
relations of the type (PQ) are clearly of countable number, we can start a contradiction argument 
similar to the abovementioned one. Therefore, there exists a resonance relation of the type 
(PQ) and denoted here by (RR), a domain Q with C e+1 boundary, £ > 3 and e > 0, such that, 
for every u G W e+2,00 (Q, M. 3 ) with ||w| 1^+2,00 < e, the domain Q + u verifies (RR). Moreover, 
since Theorem 11.21 holds true, one can assume that the eigenvalues involved in (RR) are all 
simple for Q + u with ||u||^+2 l00 < e. We then take the shape derivative of (RR) but we 
are unable to derive any contradiction. Assuming thus that this shape derivative is equal to 
zero for Q + u with ||m||^+2,oo small enough, we again differentiate the shape derivative of 
(RR) at u — 0. We then consider the variations u £ ^ Xe introduced previously and embark into 
the characterization of the main term of the second shape derivative of (RR). After lengthy 
computations (where an extra shape derivative is performed and this justifies the extra degree 
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of regularity of the boundary as compared with the argument for generic simplicity), we get 
a contradiction and conclude. It is interesting to notice the following difference between the 
proofs of Theorem 11.21 and Theorem 11.11 respectively. Indeed, for the first result, one uses, 
in the contradiction argument, the parameter defined by the angular part between x and x e 
whereas for the second result, it is the radial part between x and x £ which plays a crucial 
role. Both parameters actually result from the vectorial character of our variations and that 
enables one to adequately address the fact that (SDq) is a system of PDEs. Therefore, one 
should emphasize the flexibility of the approach proposed in this paper, which can be applied 
to genericity questions for other systems of PDEs. 

Before passing to the plan of the paper, we would like to make two remarks. The first 
one regards the reference [2D], which provides the best update for genericity questions related 
to PDEs, where genericity is meant with respect to the domain Q. Moreover, many new 
genericity results are proven there and in several situations, the author (essentially) arrives 
to the same critical issue as the one explained previously, i.e., equations of the type ([3]). D. 
Henry's approach is based on transversality theory (see for instance Theorem 5.4 in [2U1 P- 
63] instead of shape derivation. In his arguments, the main issue to contradict, "v4^ +1 is not 
dense in Ak" translates into the fact that a certain operator acting on the space defined by 
u ■ n (i.e., the domain variation restricted to dVt) is actually of finite rank. In the spirit of 
pseudo-differential calculus, he pursues the argument by evaluating that operator for functions 
u ■ n with rapidly oscillations of the type j(x) cos(u9(x)) where u tends to infinity (see [201 
Chap. 8]. The asymptotic analysis is therefore performed in terms of the phase u. Even 
though the contradiction arguments follow two different view points, our approach and that 
of Dan Henry both consider parameterized appropriate domain variations and an asymptotic 
analysis with respect to the parameters. One must however notice that the technical details 
to be handled in D. Henry's approach are much more elaborated and complicated compared 
to ours. 

We close this introduction with a brief scope on the present work and, more generally, 
on genericity problems regarding differential operators (not necessarily self-adjoint) admitting 
only a discrete spectrum. The fact that a property is generic with respect to the domain 
translates geometrically as a transversality property but the technics presented in this paper 
only rely on real and functional analysis tools. We do believe though that they are flexible 
enough to tackle similar problems for other systems of PDEs. For simplicity, we assume that 
the differential operators are completed with Dirichlet boundary conditions. Assume that 
the following two facts hold true. Firstly, one has a result of unique continuation for the 
associated overdetermined eigenvalue problem, i.e., the zero function is the unique solution of 
the eigenvalue problem where both Dirichlet and Neumann boundary conditions are satisfied 
over the whole boundary. Secondly, there exists a "good" representation formula for the 
Dirichlet-to- Neumann operator associated to the original operator, i.e., one is able to compute 
with "good" precision the kernels of the required layer potentials. Since we only perform 
computations localized in an arbitrarily small neighborhood of any point of the boundary, 
one should (in principle) only need Taylor expansions of the appropriate kernels rather than 
their explicit global expression. Once these two sets of information are available (unique 
continuation result and "good" representation formula), our technics can come into play for 
genericity questions. 

The paper is organized as follows. In Section 2, we present the necessary material on the 
Stokes system, shape differentiation and the result displayed in Eq. (J2j) and first established in 
[2"5] . The third section is devoted to the proof of Theorem 11.21 assuming that the expansion of 
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a Dirichlet-to-Neuman operator in terms of inverse powers of e is available. Then, in Section 4, 
the argument to achieve such an expansion is provided using technical results on representation 
formulas for Dirichlet-to-Neuman operators gathered in Section 6. The proof of Theorem 11.11 
is given in Section 5. Background materials on layer potentials and integral representation 
formulas for the Stokes system as well as the proofs of computational lemmas are gathered in 
Appendices A and B. 

Acknowledgements The authors would like to thank E. Zuazua and J. C. Saut for having 
suggested the problems. 

2 Definitions and preliminary results 

We start by defining precisely in Section I2~T1 the topology for the set of domains in Mr 1 with C 
boundary, where d, £ > 2. The material is standard and borrowed from (20] and [29]. We then 
recall in Section 12.21 the definition of the Stokes operator and its spectrum. The presentation 
adopted in this section is inspired by [H)J Chapter II], [3TJ Chapter 5] and [25] • Results on 
the regularity of the eigenvalues and eigenf unctions of the Stokes operator with respect to 
domain variations are derived in Section 1231 and essentially based on [22~| Chapter 7] and [2]. 
Section [23] is devoted to the shape differentiation for the Stokes system following the strategy 
of [29]. We finally recall in Section l2~5l J. H. Ortega and E. Zuazua's result obtained in [25] and 
provide an alternative proof. This result will be the starting point of our proof for Theorems 
PI and PI 

2.1 Topology on the domains 

In this section, we provide the basic definitions needed in the paper. We work in this section 
in Mr, d > 2, even though we will only be interested by the case d = 3. A domain Q of M. d , 
d > 2, is an open bounded subset of M. d . We provide now the standard topology for domains 
with a regular boundary. For £ > 2, the set of domains Q of M. d with C e boundary will be 
denoted by 3 d . Following [29], we can define a topology on B) d . Consider the Banach space 
W +1,oc (f2, M. d ) equipped with its standard norm defined by 

IMh+i,oo := supess{||-D Q u(x)||; < a < I + 1, x G O}. 

For VI eB d ,u G W e+1 >°°(n,m. d ), let VI + u : = (Id + u)(fi) be the subset of points y G R d 
such that y = x + u(x) for some x G fl and dQ + u := (Id + u)(dfl) its boundary. For e > 0, 
let V(£l,e) be the set of all f2 + u with u G W l+l,co {VL,W i ) and ||«|| w «+i,<x. < e, small enough 
so that Id + u : Q — > (Id + u) (Q) is a diffeomorphism. The topology of is defined by taking 
the sets V(Q, e) with e small enough as a base of open neighborhoods of Vt. 

A. M. Michelleti in [21] (and also reported in [201 Appendix 2]) considered a Courant-type 
metric, denoted di + \ in this paper, so that V(fl,e) is metrizable and each (V(0,,e),de + i) is 
complete and separable. For any domain Q G Of, we use Of(fi) to denote the the set of 
images (Id + u)(Q) by u G W i+1 '°°(9,, M d ), which are diffeomorphic to 11 Then is a 

Banach manifold modeled on u G W £+1 '°°(dQ,M d ) as proved in [2DI Theorem A.10]. In the 
sequel, we will sometimes identify, without further notice, the neighborhoods V(Q, e) with the 
corresponding open balls of W e+1 '°°(Q, W 1 ) centered at 0. 
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Definition 2.1. We say that a property (P) is generic in 3 d if, for every Q G H)f, the set 
of domains of Df(f2) on which Property (P) holds true is residual i.e., contains a countable 
intersection of open and dense subsets of Df(fi). 



2.2 Spectrum of the Stokes operator with Dirichlet boundary con- 
ditions 

The presentation here is inspired by [T5J Chapter II], [3TJ Chapter 5] and [23]. Let Q be a 
domain of Mr, d > 1 with C 1 boundary. We use and V(Q) to denote respectively the 

space of C°° functions with compact support in Q and the space of distributions on Q. The 
duality bracket will be denoted by (-, -)x>'xx>- 

Consider the following fundamental functional spaces for the Stokes system: 

V(Q) := {v G (H^(n)) d | div v = 0}, 

H(Q) := {v G (L 2 (n)) d | div v = in ft, u ■ n = on <9fi}. 
The space V(fi) is equipped with the scalar product of (^(fi))^ defined by 

<«,.>v:=/V«.V,:=5:/ s - s -<b, (4) 



for m := (w 1 , • • • and t> := (f 1 , • • • , v d ) in The space H{VL) is equipped with the 

scalar product of (L 2 (Q)) d which will be denoted by (•, •)#. Note that V(Q) and H(Q) are 
separable Hilbert spaces as they are closed sub-spaces of respectively (i^fi))* 2 and (L 2 (Q)) d . 
We use Lq^H) to denote the subspace of L 2 {Vt) made of the functions / with zero mean, i.e. 

f(x)dx = 0. 



d 



Remark 2.1. If we define V(ft) := {v G (V(ty) d | div v = 0}, one can show that V(Q) is the 
closure of V in (iP(f2)) d (cf. [33| Theorem 1.6, page 18]), and H(Q) is the closure of V(Q) in 
(L 2 (tt)) d (cf. [321 Theorem 1.4, page 15] and [111 Theorem 2.8, page 30]). 



Let f £ H. Since the linear form on V(f2) defined by £(t>) := / f • v, for t> G V(fi), is 

in 

continuous, by Lax-Milgram's Theorem, there exists a unique w G V(fi) such that, for every 
v G V(Q), (w,v)v = i(v) and \\w\\y < C(f2)||/||#, where the constant C(Q) only depends on 
Q. Therefore, the linear operator L from H(Q) to H(Q) defined by Lf = w is continuous. 
As L is also self-adjoint and compact (cf. Theorem IX. 16, page 169]), then, by classical 
spectral theory (cf. [U Theorem VI. 11, page 97]), the operator L admits a non- increasing 
sequence of positive eigenvalues (/ii)« e N tending to 0, and the corresponding eigenfunctions 
(^i)ieN can be taken so that they constitute an orthonormal basis of H. In particular, one has 



V0i • Vv = Xi (j)i-v, V v G V, (5) 
n Jn 

where A, := — . Note that (Aj) i6 N is a non-decreasing sequence tending to infinity. We use 

m(A) to denote the multiplicity of the eigenvalue A. 
For v G V, Eq. (j^J) is equivalent to 

(A0j + Ai0j, u )t>'xv = 0. (6) 



S 



Theorem 2.1 (de Rham-Lions). Let q G (V'(Q)) d such that 

(q,v)wxv = 0, VveV. (7) 



Then, there exists p G T>'(Q) such that q = Vp. As a consequence of Theorem \2.1\ one deduces 
from Eq. that there exists pi G V'{VL) such that 

A& + \4i = Vpi. (8) 



Remark 2.2. Note that p in Theorem 12. II is unique up to an additive constant. 

Remark 2.3. Theorem 12.11 is a consequence of a more general result due to de Rham (cf. [T0| 
Theorem 17', page 95 ]). The version adopted in Theorem 12.11 is due to Lions, also stated in 
[33l Proposition 1.1, page 14]. A constructive proof can be found in [30] . 

Remark 2.4. There exists an equivalent presentation of the eigenvalue problem for the Stokes 
system based on the Stokes operator Tg, which is defined as the operator defined on Vr\W 2 (Q) 
by Tgu G H being the unique element satisfying 

Au + Tgu = Vp, 

for some harmonic pressure field p, cf. [13 Chapter II]. Then, one has T$ = —VA where V is 
the Leray projector. One then proceeds by standard functional analysis arguments. 

The following regularity result holds for 0j and pi (cf. [33| Section 2.6, page 38]). 

Theorem 2.2 (Regularity). // the domain Q is of class C , for an integer i > 2, then, for 
i G N, (j>i_e H £ (U) and Pl G H l - l {Vt). If Vt is of class C°° , then, for i eN, foe C°°(Q) and 

We now summarize some computational results related to the Stokes system. We start 
by providing several notions of "normal" derivatives used in this context. If <fi = (<fr l )i<i<d, 

the Jacobian matrix of d> defined as (— — )i<j 7 <d will be denoted by VS. We use n to denote 

dxj ~' J - 

the outward unit normal to dQ and the superscript T used below denotes the transpose of a 
matrix. The corresponding normal derivative is given by 

^:=W. n , (9) 

and we also have 

H := (V0 + V0 r ) • n. (10) 
Finally, the conormal derivative |^ on dVt is defined as follows 

d(j) .90 

Moreover, we will use n x or n(x), with x G dQ, to denote the value of the outward normal 
vector at point x. 
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Definition 2.2. For a and b are C 1 functions defined on an open neighborhood of Q, we use 
Va : V6 to denote the following function 

Va : Vb = ^(Va + V T a) ■ (Vb + V T b), 

where • is defined in Eq. (j4j) as the Hadamard product of two matrices. 

We recall the following Green's formulas (cf. [231 page 53]). 
Lemma 2.3. Assume that d = 3. The formulas 

db 



nil 



dv 



Va:Vb + / a - (Ab - Vg), 



(12) 



and 



da 



a ■ ((A + rf)b - Vg) - / 6 ■ ((A + r/)a - Vp) 



(13) 



a- — - 

'an 0f ion ft/ j n j n 
hold for every rj £ R and for every pairs (a,p) and (b,q) of C 1 functions defined on an open 
neighborhood ofQ, taking values m R 3 x R and satisfying div a = div 6 = 0. 

Remark 2.5. One has noticed that the dot "•" has been used for scalar product a well as for 
the Hadamard product in Eq. (j4j). We will make that abuse of notation throughout the paper. 

We also need the following obvious result. 

Lemma 2.4. Let d > 2 be an integer, a e (C 1 (Q)) d H (-ffg(fi)) d and Q C R d be an open 
domain of class C 1 . Then, 

_ da T ^ 
Va = — n , on oil. 
on 



(14) 



2.3 Regularity of the eigenvalues and eigenf unctions with respect 
to the shape perturbation parameter 

In this section, £ — 1 > d > 2. Let Q be a domain in 3f. We consider perturbations u in the 
space W e+1 '°° (M. d , M. d ) with its standard norm || • ||^+i i0O - To study perturbations of eigenvalues, 
we adopt the strategy described in [2H Chapter 7, Section 6.5, pages 423-425] and also follow 
the developments of [2 Section 4, pages 1541-1548]. 

Recall that the eigenvalue problem associated to the Stokes system on Q with Dirichlet 
boundary condition is given by 



(SD S 



-A0 + Vp 
div0 



p(x)dx 



X(/> 



0. 



in Q, 
in Q, 
on dQ, 



Consider any smooth map t — > T% defined for t small enough so that To = Id and T t is a 
diffeomorphism from Q onto its image Q t '■= T t (Q). Let ((f>t,Pt, At) be the solution of 



+ Vp t 
div 6t 



\4>t 



0. 



in il t , 
in tt t , 
on dfl t , 
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By Theorem E2 <j>t e {H £ (tt t )) d n (H^(tt t )) d and p t G iZ* _1 (fit) fl Lg(fi t ). 

We next turn to the variational formulation of the above eigenvalue problem. 
For every (w,q) G (HQ(Q t )) d x Lq(^), it comes 



V0 4 : Vw dj/* 



p t div(iu) dy l + / Tr(V0t)g dy l 



We set </>' := & o T t G {H e {Q)) d n (# 1 ( n )) d and p* := pt ° T t G ff^fi) n Lg(f2). Define 
the change of variables := T 4 (y) and set := w(y t ) and r(?/) := q(y t ). Then, one shows 
that (</>*, p*) satisfies the following identity 



A{t)V<f>* :Vz- p*Tr( J B(t)V2)7(t) + / Tr(5(t)V0*)r 7 (t) = / X t ftz^t), (15) 



where j{t) = det (DT t ), A(t) = ^(t) (DTf 1 )*(DTf 1 ) and B(t) = (DTf 1 )*. For fixed t (small 
enough), all these functions defined on Q are of class C £ . 
It follows that (0*,p*) satisfies 



-div(A(t)V0') + div(p*7(t)5(t)*; 

Tr(S(t)V0*; 



p\y)i{t)dy 



\ t 0* 7 (t) in ft, 

in ft, 

on dft, 

0. 



(16) 



Let Lj, be the Hilbert space equipped with the scalar product 



(0,<0) Tt = / 0(x)V'(x)7(t) dx, 



(17) 



u G L 2 : J v(x)j(t)dx = 1. We consider C(£) and V{t) the two opera- 
tors on (Hl{Q)) d given by 



and define Lq T : 



C(t)u 



7(*) 



div(A(t)Vv) 



(18) 



and 



= -Tr ((DT^yVv) . (19) 

For fixed t (small enough), the operators C(t) and T>(t) have respectively coefficients of class 
C l ~ x and of class C . The following result holds true. 

Theorem 2.5. (cf. |H Lemma 4-2]) 

1. The operator C(t) is self-adjoint with respect to (•, -)T t and C(t)~ 1 is coercive, i.e., there 
exists C > such that, for every g G //" _1 (ft), one has (g,C(t)~ 1 g) > C || g 

2. The range ofD(t) is closed and the adjoint T>(t)* ofT>(t) with respect to (-, -)y t gwen 



2T 9 (t) 



7(*) 



div(g7(*)). 



(20) 



Moreover, the null space of T>(t) is made of constant functions on ft and its range is 
equal to L\ T (ft) . 
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Using the operators C(t) and T>(t), we rewrite System ffTB"]) as 

C(t)0* + V(t)*p t = X t (j)\ in n, 

X>(t)0* = in SI, (21) 
0* = on 

Since the operator C(t) : (ifg(f2)) d — >• (iJ _1 (f2)) d is an isomorphism, we can write 

</>* + Cft)" 1 ©* (ty = AtC(t)-V*, (22) 

and since X>(t)0* = 0, one has 

v{t)c{ty 1 v{t)*p t = x t v{t)c{ty 



Thanks to the coercivity of C(t) 1 , one concludes that V(t)C{t) 1 P(t)* is continuous and 
one-to-one in the space orthogonal to the null space of T>(t)*. It follows that 

y = A t (©(^(tr^T 1 wc(*ry. 

Finally, reporting this expression of into fl22l) . we derive that 

c(t)0* + x t v(ty {vityity'vityy^ityityy 1 = x t y, (23) 

or equivalently 

= x t A(ty\ (24) 

where we have set 

A(t) := \ld-V(ty (P(t)C(t)" 1 P(t)*)" 1 P(t)C(t)- 1 

Assume now that t h-> T t is analytic in a neighborhood of t = 0. We are therefore dealing with 
the analytic perturbation problem described in [22| Equations (6.42) page 424 and (6.47) page 
426]. Indeed, (121]) shows that the t-dependent operators A(-) and C(-) are defined on a fixed 
(i.e., t-independent) Hilbert space. We also have that 

• the operator A(t) is a closed operator with coefficients of class C 1 ^ 1 . 

• The operators t i— > C(t) and t i— >■ C(t) _1 are analytic in a neighborhood of t = 0, This 
shows that the mapping t i— >■ ^4(t) is analytic in a neighborhood of t = 0. Furthermore, 
-4(t) is bounded when t is sufficiently small. 

We next prove the following theorem. 

Theorem 2.6. Let Q C M. d be an open bounded domain of class C . Assume that X is an 
eigenvalue of multiplicity m(A) = h of the Stokes system with Dirichlet boundary condition on 
the domain Q. Then, there exist h real-valued continuous functions, u >-)■ Aj(-u) defined in a 
neighborhood V of Q in W i+1 '°°(Q, M d ) such that the following properties hold, 

• Ai(0) = A, for i = 1, h; 

• for every open interval /cl, such that the intersection of I with the set of eigenvalues 
of (SDq) contains only X, there exists a neighborhood Vi C V such that, for every 
u G Uj, there exist exactly h eigenvalues counting with multiplicity, Xi(u), 1 < i < h, of 
(S , -D(h+ u )q) contained in I ; 
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• for every u G W e+1 '°°(fl, ~R d ) and 1 < i < h, consider the map 



J -> 



t K> (Af 



R 



x n Hl{n)) d x (u^ 1 n l 2 {Q)) 

$(«), pl(u) ) 



J C M an open interval containing 0, for 1 < i < h, 0-(n) := <pt,i(u) o (Id + tn) 
and p*(w) := VtA u ) ° + ^0? where 4>t,i{u) and Pt,i{u) are respectively eigenfunction 
and eigenpressure of (SD^ +tu ). Then, for 1 < i < h, is analytic in a neighborhood 
oft = 0. Moreover, the family (0 t l (n), . . . , <j) t ^{u)) is orthonormal in Hq(Q + tu). 

Remark 2.6. This result is actually the Stokes system's version of [261 Theorem 3]. It is 
important to insist on the fact that at t = the orthonormal family 



of eigenfunctions associated to A does in general depend on u and continuity of the eigenfunc- 
tions with respect to the shape parameter u does not hold true. Therefore, only directional 
continuity and derivability with respect to u can be achieved and this is the object of the next 
paragraph. 

Proof of Theorem \2.b\ From [221 Chapter 7, Sections 6.2 and 6.5], we deduce that (A t ,0*,p*) 
defined in fT2Tl) is analytic in a neighborhood Jo of t = 0. Moreover, if A = A(0) is an eigenvalue 
of multiplicity h, by applying a standard Lyapunov-Schmidt argument (cf. for instance [22j 
Chapter 7], [T9] or [20]), one gets Theorem (J2l| when T t = Id + tu, with u e W t+1 '°°{VL, R d ), 
except that the maps \l/j's take values in R x (T^Tq (0)) rf x Jq(^)- 

Since the boundary dQ is of class C l , we now want to get that the maps t H- ((/>*(n),p*(n))'s 
take values in (H l (Vt)f]Hl(Vl)) d x (H^ 1 C\Lq{Q)) and analytic (in a neighborhood of t = 0) in 
this space. It is sufficient to show the result for the map 1 1— > <fi\{u). To see that, consider the 
power series <fi\(u) = Ylk>o A^\u)t k which is convergent in Jo as element of Hq (fi)) d . Develop 
all terms in Eq. fT2"H) as power series. A trivial induction shows that A[ k \u) e HQ(Q)) d verifies 
the equation C(0)A\ k \u) = \A(0)A[ k \u)+Rk where Rq = and Rk involves the terms A^\u), 
< J < fc — 1. By a standard argument (cf. [T7l Theorem 9.19, page 243] and [TU Theorem 
5, page 323]), one shows by induction that firstly Af\u) G (H e (Q)) d (since the operators 
A(t) and C(t) of class C^" 1 ) and secondly their (H e (Q)) d - norms verify the appropriate upper 
bounds insuring that the map t £ fc>0 Af\u)t k takes values in (H e (n)) d and is a convergent 
power series in some neighborhood J of t = contained in J . 



2.4 Shape differentiation 

The subsequent developments follow a standard strategy (cf. [29], Theorem 2.13] for in- 
stance) but seem to be new for the Stokes system with Dirichlet boundary conditions. Fix 
u G W e+1,oc (Q, R ) and set T t = Id + tu for t small enough. In this section, we define and 
calculate the differential systems verified by the derivatives at t — of the eigenfunctions 
((f)i } t(u) , Pi,t(u)) defined in Theorem 12.61 For that purpose, we must first consider the deriva- 
tives of the maps <p\{u) and p\(u). Since we perform such a computation along a fixed analytic 
branch (\\(u),(p t> i{u),p t ^{u)), the index % is omitted for the rest of the paragraph. 



(0o,i( M )> • • ■ >0o,h(w)), 



□ 



13 



According to Theorem I2.6[ is analytic in a neighborhood of t = and we set 

cZ0*(it) dp* (it) 



i=0 



«„) := ^ . , (25) 



t=0 



We next proceed in a similar way as in [291 Theorem 2.13]. For every open set oj whose 
closure is included in Q, we consider (0t(u))| w and the restrictions of 4>t{u) and pt(w) 

respectively to oj. As compositions of two analytic maps in a neighborhood of t — 0, (&(«)) | w 
and ^(m))!^ are also analytic in a neighborhood of t = and their derivatives at t — are 
equal to — V0 • and (p(u) — Vp ■ respectively. It is then easy to see that these 

formulas are actually valid over the whole Q and thus, if we use <fi'{u) and p'(u) to denote the 
derivatives at t = of <pt and pt respectively, one finally gets that 

4>'(u) = 4>{u) — V0 • u, p'{u) = p{u) — Vp ■ m, in Vt. (26) 

We refer to 4>'(u) and p'(u) as the shape derivatives in the direction u of the eigenfunction and 
eigenpressure (4>,p) associated to A. 

According to Theorem ESI <fi(u) and p(u) belong to {H\Vt) n H^(Q)) d and H l ~ x n Lg(fi) 
respectively. Therefore, they admit traces on <9f2 in H e ~ 1 ^ 2 (dVl) and H e ~ 3 ^ 2 (dVl) respectively. 
Since i — 3/2 > d/2, these traces are continous functions on dfl, which are of course equal to 
zero. From Eqs. (126]) and ffT6]) . we deduce at once, by using Eq. ffT4l) that p'( u ) + div(wp) e 
Lg(ft) and 

4>'(u) + (w • n)— = on dQ. 
on 

It remains to determine the relations satisfied by the derivatives <f)'(u) and p'(u) inside the 
domain Q. For that end (see [21 Proposition 4.6] for more details), we take the derivative 
with respect to time evaluated at t — of Eq. . For arbitrary test functions (z, r) e 
(V(Q)) d x V(Q), we obtain 

J (A'(0) V0 + V0(m)) : Vz - J (p{u)div{z) + pTr{B'(0) Vz) + pdiv(z)y (0) 



+ / (Tr(5 / (O)V0) + div(0( M )) + div(0) 7 '(O)jr = / (V(w)0 + A0(w) + A0 7 '(O) J 2. (27) 

To simplify the previous equation, we use the following relations between time derivatives and 
shape derivatives, 

V(0) = div(u), A'(0) = div(u)Id - (Vu + V T u) and B'(0) = -V T u. 

We first use the boundary conditions for <fi and notice that the term multiplied by t'(O) in the 
integrand of Eq. (j2"7|) is the PDE satisfied by <fi. Eq. 027|) then reduces to 



J (^V0'(m) + V(V0-m)-(Vm + V t m)V0) : V2~y (( P '(u) + Vp-u)div(z)-pTr(V T uV 

\'{u)(j) + \(j)'{u) + AV0 • uj z, 

y ( - Tr(V T nV0) + div(0'(u)) + div(V0 ■ u))r = 0. 



and 
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After some integrations by parts and using the boundary conditions, one deduces the two 
identities 

f V0'(m) :Vz+ I Vp'(u)-z = [ (\'(u)4> + \(f)'(uj)z and f div(0'(«))r = 0. 
Jn Jn Jn ^ ' Jn 

These identities hold for every (z, r) G (T>(Q)) d x T>(Q), and they yield to the equations which 
are valid in Q 

-(A + \)<p\u) + Vp'(u) = -A'(u)0, div(0'(u)) = 0. 

In summary, the shape derivatives <fi'(u) and p'iu) satisfy the following inhomogeneous Stokes 
system of PDEs 



f -(A + A)0'(tO + Vj/(u) 
div 

0(u) + (wn)^ 
p'(u) + div(-up) 

2.5 Ortega-Zuazua's result 



= -A'(m)< 




G Lg(fi). 



in f2, 
in f2, 

on 



(28) 



Our argument for establishing Theorem 1 1 . 2 1 requires the shape differentiation of the eigenvalue 
problem (SDq). The first step of the contradiction argument (i.e., assuming that the simplicity 
of the spectrum is not generic) was already conducted by J. H. Ortega and E. Zuazua in |25j . 
We next recall precisely the main result they obtained and, for that purpose, we introduce the 
following definition. 

Definition 2.3. Let £, d be two integers such that £ — 1 > d > 2. A domain Q G verifies 
Property (Poz)d if; f° r every A eigenvalue of the Stokes operator with Dirichlet boundary 
conditions (SDq), one has m(A) < d — 1 and if m(A) = d — 1, for 1 < i, j < d — 1 and i ^ j, 
the following three conditions must hold on dQ, 



On 



ii 



dn dn 



0. 
0. 



d(f)i 




<% 


dn 




dn 



(29) 
(30) 
(31) 



where the 0j's, 1 < % < d — 1 are orthonormal eigenfunctions associated with A. 
Then, the main result in [23] is the following. 



Theorem 2.7. Let £,d be two integers such that i — 1 > d and d is equal to 2 or 3. Then 
Property (Poz)d defined above holds true, generically with respect to Q G V>f . 

As an immediate corollary, it is proved in [25] that Property (Simple) holds true generi- 
cally for domains in M. 2 . Since we adopt a viewpoint different from |2"5] , we provide below a 
complete argument. We need to provide the following definition, similar to that of "minimal 
multiplicity" in [201 P a g e 56]. 
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Definition 2.4. Let Q G Bj and A an eigenvalue of (SDq). We use mn(A) to denote the 
liminf over the multiplicities m(A n ), where A n is an eigenvalue of (SDn n ) such that Q n — > Q 
and A n — > A as n tends to infinity. 

Several remarks are now in order with regard to Definition 12.41 

Remark 2.7. There exists a sequence of domains (f2 n ) in B d e and a sequence (A n ), where A n is 
an eigenvalue of (S^i n ), such that fl n — > fl, A n — > A as n tends to infinity and m(A n ) = mn(A) 
(and it is also equal to mn n (X n )). 

Remark 2.8. Moreover, Property (Poz)d for a domain Q G 3 d is clearly equivalent to the fact 
that, for every A eigenvalue of (SDq), mn(X) < d — 1, with the equality case described by 



and, for I > 1, consider 

^ := {f] +u 64tie W^ +1,oc (n ,R d ), m^ (A) < cZ-1 for the first I first eigenvalues of [SD Qo+u )}. 
Set A := f] leN Ai. Note that 

A = {Q +u e Aq, ue W e+1 '°°(tt 0j R d ),m Qo+u (X) <d-l if A is an eigenvalue of (SD Qo+u )}. 

The proof is based on the application of Baire's lemma to the sequence {A;}^^. As Ai is 
open in A for every I G N, we only need to prove that, for / G N, A+i is dense in A\. 

We proceed by contradiction. Assume that A[ + i is not dense in A\. Then, there exists 
u G Ai \ Ai + i and a neighborhood U of u such that U C Ai \ Ai + ±. Set Q := Q + u and let 
A be the (7 + l)-th eigenvalue of (SD^). For s > 1, let A s (-) be the function which associates 
to 6 the s-th eigenvalue of (SDq). Note that A s (-) is continuous and A = A;+i(f2). 
According to the contradiction assumption, one has m := m^(X) > d and then Xi(Cl) < X. 
As a consequence, if (f2 n ) is the sequence in considered in Remark 12.71 and associated to 
Q, then it has the following additional property: for n large enough, there exists e n > such 
that, for every Q' with d(Q', Q n ) < e n , one has that 



In particular, m(A; + i(-)) is locally constant, equal to m > d in an open neighborhood of Q n , 
for n large enough. We will contradict that latter fact, i.e. the existence of a domain fi* 
where m(A; + i(-)) is constant and equal to m > d in an open neighborhood U* of For 
simplicity, A is used to denote A^+i^*) in the remaining part of the argument. Once for all, 
fix an orthonormal family v = (v±, . . . ,v m ) of eigenfunctions of (SDn t ) associated to A and 
define the m x m matrix 



Note that M(v ) is real symmetric. We next perform shape differentiation with respect to the 
parameter u G [/*. Using the notations of Theorem I2.6[ we consider, for every u G [/*, the m 



Eqs. d22D,(I3DD,(ISI]). 



Proof of Theorem\2. 7[ Fix a domain Qq G D^. We define, for / G N, the sets 



A :=B d e (Q ), 



m(A z+ i(n')) 



m > d. 
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analytic branches t (A*(u), <p t ,i{u) , p t ^(u)) , for % — 1, . . . , rn, given by Theorem 12.61 We use 
(j){u) := ((j)i(u), . . . , 4> m (u)) and (pi(u), . . . ,p m (u)) respectively to denote 

(00,1 0), • • • , 4>0,m(u)), Oo.lO), • • • , 90,m(w)), 

the eigenfunctions and eigenpressures associated to A (i.e., which correspond to the values of 
the t) i(«)'s and ^(i^'s at t = 0). 

Since v and <p(u) are orthonormal families of eigenfunctions associated to the same eigen- 
value A, then, for every 1 < i < m, there exists S(u) £ SO{m) such that (p{u) = vS(u) (with 
the convention that the 0j(w)'s and the v^s are viewed as column vectors of R m ). One clearly 
obtains that 

M((j)(u)) = S(u)M(v)S(u) T . (32) 
We now need the following standard result whose proof is given in Section iB.ll of Appendix. 
Lemma 2.8. Using the notations defined above, then 

diag(A:(n)) 1 <,< m = -M{<j>{u)) (33) 

holds for every u £ W e+1 '°°(n,R d ). 

We next proceed with the proof of Theorem 12.71 

The fact that m(A^ + i(-)) is constant and equal to m in a neighborhood of u = is equivalent 
to the fact that \\{u) = A*(m), 1 < i, j < m, for t small enough, implying that A'j(w) takes 
only one single value /i as i runs from 1 to m. In other words, M(0(w)) = — /ild m and then 
one gets 

M{v) = -//Id m , 
thanks to Eq. fl32l) . That yields the equations 

2 



[u ■ n 



an, 



X 



dvi 


2 




dn 




dn 



,ii ■ n 



an, 



dvi dvj 
dn dn 



0, for 1 < i,j < m, 

0, for 1 < i,j < d — 1, i 7^ j. 



(34) 
(35) 



The integrals in the above equations define linear maps in (u.n) and are equal to zero in an 
open neighborhood of u — 0. It thus implies that, for distinct 1 < i,j ' < m, 



dvi 




dvj 


dn 




dn 



on dfl* 



"7t ■ =0 on 9i2*. 

an an 

Moreover, using Lemma (12 .4p . one has, for 1 < i < m, 

— — ■ n = 0, on Sfi*. 
an 



(36) 
(37) 

(38) 



Assume now that there exists Xq £ Sfi* and an index i £ {1, • ■ • , m} such that — — (xq) is not 

an 

zero. According to Eqs. ( 1361) . ( 1371) and ( )38|) , the (m + 1) vectors given by — —(x ), 1 < j < m 

dn 
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and n(xo) are all non zero and two by two perpendicular. This is a contradiction because 

3v' 

these vectors belong to a <i-dimensional vector space. Therefore, must be identically equal 

on 

to zero, for 1 < % < m. 

Thanks to a unique continuation type of argument due to Osses (cf. [27]) and which is 
valid only for d equal to 2 or 3, one concludes that the v^s must also be identically equal to 
zero, which is in contradiction with the facts that the v^s have L 2 -norm equal to one. 

□ 

Remark 2.9. This argument is an adaptation of the original proof by J. H. Albert in [TJ to the 
Stokes system with Dirichlet boundary conditions, and the perturbation parameters being the 
domains of M. d . See also [20| Example 4.4] for a more general situation. 



3 Proof of Theorem 1.2 



For the rest of the paper, domains Q are bounded subsets of M 3 with C boundary, i.e., 
d = 3 and I > 4. We follow the classical strategy initiated by J. H. Albert in [1] for the 
Laplace operator with Dirichlet boundary conditions. This strategy was in particular applied 
successfully in [23] for the generic simplicity of the Stokes operator in two space dimensions, 
and in [9] for other Laplacian-like operators. Fix a domain Q G D 3 . We define, for / G N, the 
sets 

A := D 3 (ft ), 



and, for / > 1, 

At := {fio+w, u G W^ +1 '°°(ft ,K 3 ) 

Set A : = fl/eN^- Note that 
A = {u G W^ +1 '°°(fio,M 3 ), Q 



flo+u G A and the I first eigenvalues of (SD^ 0+U ) are simple}. 



u G A Q and the eigenvalues of (S-Dqo+u) are simple}. 



Again, the proof of the generic simplicity of (SDq) is based on the application of Baire's 
lemma to the sequence {Ai}i e ^. As Ai is open in D|(Q ) for every I G N, we only need to 
prove that, for / G N, A+i is dense in A[. We proceed by contradiction. Assume that Ai + i 
is not dense in A\. Then, there exists u G Ai \ Ai + \ and a neighborhood U of u such that 
U C Ai \ Ai + i. By Theorem 12. 7\ we can assume, without loss of generality, that there exists 
Q := Q + u for some u G U verifying the following: there exists an open neighborhood 
V dU of such that, for every u G V, then Q + u verifies: 

(i) the first I eigenvalues Xi(u), . . . , Xi(u) of (SDn +u ) are simple; 

(ii) the multiplicity of the (I + l)-th eigenvalue Xi + i(u) of (SD^ +U ) is equal to 2 and, on 
dQ + u, one has 



On,, 



dn u dn u 



0. 





i = l,2, 



901 




d(f) 2 









(39) 
(40) 
(41) 
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where n u is used to denote the outer unit normal at dVL + u and (fa, fa) is any pair of 
orthonormal eigenfunctions associated with A; + i(it). 

Remark 3.1. These conditions simply state that, for an eigenvalue A of (SDq) (say the (1 + 1)- 
th), its multiplicity is larger than or equal to 2 and, for every variation v in W +1 '°°(fi + tt, M 3 ), 
there are two equal directionnal derivatives (in the direction of v) of A/ + i at u. This fact 
actually does not depend on the dimension d > 2 of the domain Q. In dimension two, the 
above conditions immediately yield that 

d<j>i _ dfa _ Q 
dn u dn u 

for any pair of orthonormal eigenfunctions associated with A/ + i(u), and one derives at once a 
contradiction by the unique continuation result of [27], see also [23]. However, in dimension 
d = 3, conditions ( |39|) . fl4"D]) . and ( |4T|) do not immediately yield a contradiction since three 
non-zero two- by- two orthogonal vectors may exist in dimension d = 3. 



3.1 Shape derivation of Equations (1391) (1401) and (1411) 

We begin with the following preliminary result. 

Lemma 3.1. The shape derivative fa of (pi in the direction V satisfies 

where we use V n to denote the normal component V ■ n of the direction V evaluated on dQ. 

Proof of Lemma VS. 1\ From the fact that fa vanishes on dQ and satisfies div(fa) = 0, one 
knows that 

{Vfafn = 0. (43) 
Taking the shape derivative of the two sides of Eq. (143|) . one gets 

{Vfaf n + (Vfa) T n' + y n A((V0,) T n) = 0. 
(dfa T y f d^y 



Since (Vfa) = n J = n \q~) > ^ comes that 



mf n + n')n + V n ^ ((Vfafn) = 



hence 



(V^ T n = -(^,n')n-V n -^((Vfafn). 



The proof is finished once we report this expression in the definition of the co- normal derivative 

of fa. 

□ 

Proposition 3.2. If fa satisfies / fffPj) and ^Ify , then we have, for j = 1,2, 



dv ' dn I \ dv ' dn 
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Proof of Proposition \3.£\ The shape derivative of Eq. fl39|) gives 



!*\W*M\ = - f »|( ?1? )l (45) 



dn 1 dn I \ dn' dn / n dn\\ <9n ' dn 



Since — — • n — 0, it comes from Lemma (13. li that 
on 



\ on on I \ dv on I \on on 
hence we deduce that 



dn dn I \ dn dn 

/ (90- d(j)j \ _|_ / ^4>i \ y f/ 9 ^ T d<f)j \ _|_ / <9 

\ dv ' dn / \ dv ' dn / n \\dn 1 dn I \dn 3 dn 
From Eq. (H5l) . we get after identification that 



A , / a< ?i ,r f d f/ d 0j 9 4>i\\ . / d T d(j>j\ / d T , 



dv'' dn I \ dv , dn / n \dn\\ dn 1 dn I J \dn 1 1 dn / \dn 3 1 dn 



and this ends the proof of Proposition 13.21 

□ 



3.2 Special choice of V n 



d<Pj 



Let x G <9f2 such that the vectors -^—(x) and -—— (x) span the tangent space TJdTl). Let 

dn dn 

1A X be a neighborhood of x in dVL such that, for all y belonging to U x , the vectors ^(y) and 
^-(y) span T y (<9f2). For ?/ G <9f2 near x, we write the parametrized form of dQ near x as a 
graph over the tangent plane at x : if rj = P x (y — x) is the orthogonal projection of y — x onto 
the tangent plane T x (dQ) with n sufficiently small, there exists an open neighborhood T X U X 
of in T x (dQ) such that the map h x given by 

h x '■ T X U X i y IA X ( A(\\ 

7] !-)■ y = x + 7] - v x (rj)n x , 

is well-defined and is a diffeomorphism onto its image. For y near x, we have 

^0?) = ^V T KxV + 0(\r]\ 3 ), asr/^0, 

where is the symmetric matrix representing the curvature operator at x. We fix once for 
all 5 > small enough so that |?7| < 26 implies that y = x + r] — v x (rj)n x belongs to IA X . 

One has the following standard relations where all the vectors are embedded in M. 3 and (•) 
denotes here the standard scalar product in R 3 . 
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Lemma 3.3. Let n' x {rj) = ^-n x {rj) G T x (dfl). We have 



n y 



ii) (n x , y-x) = -^rj T K x rj + 0(|r/| 3 ) as r) ->■ 0, 
*) (n x ,n y ) 



(47) 



= l-^|^| 2 + 0(|r/| 3 ) as V ^0. 

Proof of Lemma Iff. 31 These equations are easily obtained by standard facts from the theory 
of surfaces (cf. [7J Chapter 10]) and are explicitly given in [201 page 146]. □ 

Remark 3.2. We note that the inverse of the Jacobian of the change of variables h~ l : y — > 
rj = h~ 1 (y) from a neighborhood of x on dQ to a neighborhood of in R 2 is equal to (n x , n y ). 

We are now ready to define V n . Let e « 5 be a positive real number. For 7] G B E (0) C 
T x (dfl), consider the points x G U x which can be written as 

x = x + rjo - v x (vo)n x , 

and set 

rjo = r (cos9 ,sm9 ) T , 

where 6*0 G S* 1 and < r < e. 
Our choice for V n will be 

V n (y) := (a £ , v Moh~\y), (48) 
where, for rj G M 2 (identified with T x (dQ)), 



a s,vo(v) := ^exp 



1 r I V ~ Vo 



12 



e 2 e 2 



and is a smooth cut-off function equal to 1 on -8(0, 35/2) and on R 2 \ -8(0, 26). 

Lemma 3.4. If y = h x (rj) with t] G B(0,S), then we have 

W n (y) = W(h x {y)) - (W{h x (y)), v' x (h x {y)))n x , (49) 

where W n (y) is used to denote the tangential gradient ofV n along dQ and W(h x (y)) G T x (dQ) 
is given by 

tnMz/)) = Va £tVo (h x {y)) - — — v x {h x {y)). (50) 

Proof of Lemma \3.4\ Recall that 77 = h x (y). By definition of V n , one has, for every tangent 
vector w G T„(dfi) with f = dh~ 1 (y)w G T x (<9fi) 

(W„(y),«;) = dV n (y)w = d(a £ ^ vo o h~ 1 )(y)w = da £tV0 (7])dh^ 1 (y)w 
= da £!V0 (r])£ = (Va E , V0 (ri),£), 
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where da etVo denotes the differential of the scalar function a Em . From Eq. (|46[) . one gets that 
w = £ — (v' x (v)yO n x- By a simple computation, one deduces that 

f = w + (v'xiv),™)™*, 

where the latter inner product is taken in R 3 . According to the orthogonal sum IR 3 = T x (dQ) © 
M.n x , one has VV^(?/) = W(r)) + (V n (y), n x )n x for some vector W{rj) G T x (dfl). Since VV^(y) G 
Ty^Q), on gets after taking the inner product of the decomposition of V n (y) with n y that 

(V n (y),n x ) = -(W( V ),v' x (r ] ))- 
Plugging the two previous displayed equations into Eq. floTl) . one deduces that 

W{ri) + (W( V ),u' x ( V )H(v) = Va^fo), 
which in turns yields Eq. fl50l) and hence Eq. (j49|) . □ 

Convention For the ease of notation, the gradient of a scalar function will be considered in 
the following as a line vector instead of a column vector. This convention will allow us to use 
the notation V for both scalar and vector-valued functions in a consistent way. 

3.3 End of the proof of Theorem 11.21 

The main technical result of the paper is summarized in the following proposition. The proof 
is provided in Section HJ 

Proposition 3.5. Let x G dQ such that the vectors —^(x) and -—^-(x) span the tangent space 

on on 

T x {dVt). Use P x to denote the orthogonal projection onto T x (dVt). Then, for e small enough, 
one has, for every tjq G B e C T x (dQ) and the corresponding variation V n defined in Eq. 
that, for j = 1,2, 

P'^-W) = ^[M^(f ) + M 5 A ^f )-flM^(f ))^(x) 



+2^M^(f ){f)o, + (52) 

where fo = fjo = jr^y =: (cos(#o)j si n (#o)) T ^nd M Al (-), 2 < k < 5, are nonzero entire 
function defined in Eqs. [8~4\) , U86\) and U85\) respectively. 



We can now conclude the proof of Theorem II .21 By conditions (1291) and f l30|) . and Propo- 
sition I3.5[ we have 

Ax) ^(x)) + Ax) ^(x)) 
{ du [ h dn [ )} + { du [ h dn { >' 



ne 6 on on e 



-M 4 Al (f )^cos(^ - O ) cos(£ 2 - 9 ) + O(-), 
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with —2-(x) = r J cos 6 j, sin 8 A , for j = 1,2 
On 

However, Proposition 13.21 implies that 



f x) ?4l( x )) + i d Jk {x) ?h (x) ) = o(-) 



Therefore, if we now fix ro < 1 such that M/ 1 ^) ^ and recall that r<j> > 0, we have, for 
every 9 G S 1 , 

cos(9 1 -9 )cos(9 2 -9 ) = 0(e). (53) 
By letting e tend to zero, we deduce that 



cos(#i - O ) cos(0 2 - O ) = 0, 



(54) 



since 9 does not depend on e. Again, by conditions (|2"9"1) and (1311 . one has |#i — 2 | — 7r /2- 
Then, by replacing the arbitrary angle 9$ by 0q — 0i in Eq. (1331) . one derives that 



sin 20 o = 0, 



(55) 



holding for an arbitrary angle 9q G This yields the final contradiction and Theorem 11.21 is 
established. 



4 Proof of Proposition 13.5 



This section is devoted to the proof of Proposition 13.51 The argument starts by applying ( I162[) 
to (f) x — 0'-, j — 1, 2, solution of (I165]) - (|168|) . The four terms of the right-hand side of (I162[) 



correspond to four terms W/, 1 < i < 4 respectively. Since </>'■ 



-r, 



that 



an 



dv 



[x) = W((x) + Wi(x) + Wi(x) + (x) 



on <9f2, it comes 



(56) 



where we have in coordinates, for 1 < s < 3, and 0j = (</>™)i< m <3, 

<9 2 r° m (* - ?/) 



(57) 



E [(-2)^ 

fc=l 
v 

E [(-2)(^' 



1 fc 



k=l 



Wi(x) 



E [(- 2 )(* 



W?)(ar), 



(58) 



fc=0 



9Q dN(x)dN(y) dn 
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and 



N 



Wi{x) 



R 



£ (- 2 )(* 



p.v. 



R 



k=l 
N 

£ 

k=l 



d 2 T° (x-v) dcbf 
dn dN(x)dN(y) dn 



d 2 Aj n (x - y) 
dN(x)dN(y) 



v n(v)^(y) My)- 



(60) 



1 I h~Hy) - h~ l (x ) I 2 
We take V n (y) = — exp[ - ^ ] and tackle the asymptotic expansion of each 

term appearing in the right hand side of the equation quoted above. Our strategy is simple: 
we show that the main term of the expansion is contained in W(, where appears the effect 
of the hyper- singular operator. Next, we prove that all other terms W?(x), i = 2,3,4 are 
actually remainder terms. These are the contents of Proposition 14.11 and Proposition 14.111 
respectively given in the next subsections. 

4.1 Expansion of W{ 

The goal of this subsection is to provide the main term in the expansion of W( (x) defined in 
Eq. ( 157)) . More precisely, we prove the following. 



Proposition 4.1. With the notations of Proposition 1X5). we have, fore > small enough 
and j = 1,2, 



Px(W((x)) 



+ 2- 



-r 2 

e o 

-r 2 



M 2 Al (f ) + M 5 Al (f ) - f 2 M^(f 
M^(r )fo,^(s))fc + O(i). 



dn 



[x 



(61) 



4.1.1 Computational lemmas 

We begin by studying the term W((x) defined in Eq. (157)) . We start with the following lemma 
whose proof is deferred in Appendix. For u = (w m )i< m <3 : dQ >->■ R 3 , we will use E(u)(x) to 
denote the value at x G dfl of the hypersingular operator 



E(u)(x) 



Jan dN(x)dN{y) 



1< s < 3. 



(62) 



Lemma 4.2. Let a : dQ i— > R and ip : dQ i— y R 3 be C l functions. One has 



iirE(ai>)(x) = Ai(a,ij})(x) 



(63) 



i=l 
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where 



(n x ,n y ) 



Ai(a,^)(x) = p.v. / ^^{^{y),x-y)V 1 a{y) + {Va{y){x-y))^{y))da y , (64) 
Jdn \ x — y\ v 7 

A 2 (a,^)(x) = p.v. / a( f )K, ,3 y) (W(y) + V T ^))(x - y)^, (65) 
Jan F ~~ 2/1 v J 

a t i\i \ f ( n x,^(y))^oi(y)(x-y)-(ijj(y),x-y)Va(y)n x 

A 3 (a,i/j){x) = p.v. / j fg n^dcTy, (66) 

Jdn \ x — y\ 

At \\( \ f a (y)( n x, (V^(y) - V T ijj{y)){x - y)) 

A±{a,ip)(x) = p.v. / : - n y day, (67) 

Jdn \ x ~y\ 

A 5 {a,iP){x) = I l{x,y)[V{a^){y)]da y , (68) 
Jan 

where /(•,■) «s a weakly singular operator of class C% (l) (see Appendix \A.S\ for a definition). 



Lemma 14.21 will be applied with a = V n and ^ = — — , j = 1,2. We will consider the 

cm 

change of variables introduced in Section 13.21 and, using these notations, we set 

( cos 9 \ ( cos^q \ , . / cos6L \ _ ?7o 

V sin 6 1 / V sin J y ' v y sin ^ J e e 

Recall that, with the conventions of Subsection 13. 2\ one has f < 1. In the sequel, we will 

provide an asymptotic expansion for each of the Ai, 1 < % < 5, using powers in the variable -. 

e 

e~ r o 1 

We will have two types of terms, one of the type Xi (or XA and the other one of the 

_ s 2 

e 

type Yi, where is an integer and Xi, Y^ are vectors with bounded norms. For each Ai, 

e mi 

e _f o 1 

1 < i < 5, we will identify the term of the first type (i.e., Xi or XA with the largest 

value of mi, then gather them and consider all the others terms as a rest. For that purpose, 
we will use repeatedly the following two lemmas whose proofs are deferred in Sections IB.3I and 
IB. 41 in Appendix B. 

Lemma 4.3. With the notations above and for any non negative integer m, one has 

< W (69) 

1—771 ' — £ m 



'5(0,(5) I V 

with C(m) a positive constant only depending on m. 
Lemma 4.4. With the notations above, 

p . v ./ ^Mn in = tl M tm%, (70) 

Jr 2 I V I £ 
where M^ 1 ^) is a nonzero entire function defined in M8$fy or U83\) below. 
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We will provide detailed computations for Ai(a,ip)(x) in the expansion of W((x) and will 
only sketch the main steps for the other terms. In these computations, we will systematically 
refer to the following procedures. 

(PI) The first one consists of decomposing a C 1 vector-valued function F(y) in two parts as 
F(y) = F(x) + G(x)(y — x), where G is a continuous matrix- valued function. 

(P2) The second procedure consists of cutting an integral / • • • da y as 

Jan 



■ ■ ■ d<jy — I ■ ■ -drj — ■ ■ ■ drj + / ■ • • drj, 

dn J B '(0,2(5) JB(0,6) J B(Q,2S)\B(p,S) 

_ s 2 

e ^ 

and majorizing the second one by Ci for appropriate constant Ci and integer wij. 

£ m.i 

Finally, note that (ip(x),n x ) = 0. 

(P3) In certain integrals of the type J dn ---da y = J B ( 02 s) ' " " dr}, the term W n (y) will be 
expressed after changing variables as 

Va^oO?) + G(r])X7a £ , V0 (r]) + H(rj)n x , (71) 

where G(rj) denotes the non-positive symmetric matrix — ~fffi~7fa^ and H(rj) is a real- 
valued function. Note first that, for r\ small enough, one has that |G(^)| < C\rj\ 2 for 
some universal positive constant C. Thus the contribution arising from G(?7)Va £ Vo (v) i n 

dk 1 
the expression of -^-{x) will be shown below to be trivially a O(-). Moreover, in the in- 
tegrals abovementionned, only their contributions tangent to T X (Q) are relevant, thanks 
to Eqs. f l29|) and ( jjjj) . In conclusion, it is enough to only estimate the contribution of 
(rj) in Eq. ([7T]). 

4.1.2 Asymptotic expansion of A 1 

We give in this paragraph the asymptotic expansion of Ai(a,i/})(x) with respect to e. Recall 
that 

A 1 (a^)(x) = p.v. [ i!^((^(y),x-y)V T a(y) + (Va(y)(x-y))^(y))da r 
Jdn \ x — y\ v J 

Proposition 4.5. For e > small enough, one has 

-f 2 

P x {AiM){x)) = 2^(M 2 Al (f ) + M 5 Al (f )-f 2 M 3 Al (ro))^(x) 



2^M^(f )(r ]0l ij(x)}r ]0 + O(i). (72) 



For the sake of clarity, we set Ai(a,ip)(x) := Ai t i(a,if))(x) + Ai t2 (a,ip)(x) with 

A lt i(a,i/))(x) := p.v. / p^My), x -y)Vr a (y)da v , (73) 

Jdn F — V\ 

A 1>2 (a,i;)(x) := p.v. / |^(Va(y)(a: - v))1>{y)da v . (74) 

Jan F — v\ 

We will establish separately estimates of these two terms in Lemma 14.61 and Lemma 14.81 
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Lemma 4.6. For e > small enough, one has 

P x (A hl (a,4>)(x)) = 2^(M 4 Al (f )(r/ ,^(a;))r/o + M 2 Al (fo)^(x)) +0(^), (75) 

where M Al (-) and M Al (-) are non-zero entire functions defined by (Q?2j) and ( fffSj) respectively. 

Proof of Lemma \4-6\ Using the change of variables introduced in Subsection 13.21 and taking 
into account Lemma [3.41 and Remark 13. 2[ we have 

PMvMM) = \ p v. f "«'"><"- ^">"-^M> (i d2 + G( ,)) ( , f - »)*,. 

Then, by taking into account Procedures (PI) and (P3), 

P x {AiM^){x)) = I A ^(a^)(x) + J A ^{a^){x) + R A ^(a,*j;), 

with 

/ Al,1 (a,^) := — p.v. / ! — — (77-770)^77, (76) 

£ JB(0,<5) M 

J^(o,«(«) := i / (77) 
P^( a ^)( x ) := /" (78) 

^B(0,2<5)\B(O,5) 

where, in P" 41 ' 1 (a, ^)(a;), one has the same integrand (in local coordinates) as in Ai(ot, V'X 2 ')- 
Clearly, there exists a positve constant C$ only depending on 5 such that, for e small enough 
with respect to 5, one has 

HP^a^X*)!! <C^. (79) 



5 



Moreover, one can apply Lemma [4.31 to J 1,1 (a, i[>)(x), one gets that 

ll^ 1 («^)(*)ll<^(C(l) + ^), 
and since — = 0(1), one finally deduces that there exists a positive constant C* such that 

\\J A ^(a^)(x)\\<^. (80) 

Note that, for e small enough the upper bound of fISDl is larger than that of (J7SJ). 
It remains to estimate I 1,1 (a, if))(x). First of all, notice that the norm of 



s 2 



K 2 \B(0,<5) \V\ 



C$e 3^ 

is clearly less than or equal to for some positive constant C$ only dependent on 5 and 

e small enough with respect to 5. 
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We can therefore estimate, instead of I AlA (a, ip)(x) : the quantity J" 41 ' 1 (a, ili){x) defined by 

I \ot,%l)){x) := — p.v. / — '- — — (r) - r) )drj. (81) 

By using polar coordinates, one gets 
I A ^{a^){x) 

e" f o r°° r 2 f 2w r_ /cos6>\ 

= 2^^r<p exp( — -)dr / cos(6» - 6^) exp(2-r cos(6> - 6» )) . „ c$ 
£ Jo £ Jo £ \smoj 

e-'o /cos0 o \- p exp(-g) ^ ^ r_ 

-2 — — r^, I . I r p.v. / — dr / cos(0 - exp(2-r cos(0 - Vo))dv 

e \smu j Jo r Jo £ 



-ry, 

where 



e" f o / M Al (f ) cos(0 o - 0^) cos # + M 2 Al (f ) sin(0 o - ^) sin O 
"e 3 "^ ^M^ 1 (f ) cos(0 o - Oj,) sin O - M^ 1 (f„) sin(0 o - ^) cos O 

-2^cos(0 o - r 2 M 3 ^(f ) 

, e" f o / [M x Al (f ) - f M 3 Ai (f )] cos(0 o - 4 ,) cos O + M 2 Al (f ) sin(0 o - 0^) sin 6 ' 



V [M^ 1 (f ) - foM 3 Al (f )] cos(0 o - 9*) sin O - M 2 Al (f ) sin(0 o - ^) cos O , 



M Al (f ) := / exp(-r 2 )dr / cos 2 0exp(2rf o cos0)d0, (82) 
Jo Jo 

1*00 f'l'K 

M Al (f ) := / exp(-r 2 )dr / sin 2 0exp(2rf o cos 8)d6, (83) 
Jo Jo 

1 f°° exr>(—r 2 ) f 2n 
M Al (f ) := —p.v./ — ^ ^dr / cos0exp(2rf o cos0)d0.. (84) 

r o Jo r Jo 



The needed information about the functions M i 1 (•),«= 1, 2, 3, is gathered in the following 
lemma, whose proof is given in Section IB.5I in appendix. 

Lemma 4.7. Fori = 1,2, M Al (-) are entire functions . Moreover, the function M Al (•) defined 
by the relation 

Mt{z) := ^(M^(z) - z 2 M^(z) - M£{z)) (85) 



1 

z 

is a nonzero entire function. 



Using Lemma [4 .7\ we further simplify I 1 11 as follows. 

I x 1,1 (a, V) 



2 e- f g r / [M^(fo) - rgM*(r )] cos(0 o - fy) cos0 o + M^(f ) sin(0 o - fy,) sin 0c/ 
e 3 H UM^(ro) - r 2 M^(f )} cos(0 o - 0yO sin0 o - M^(f ) sin(0 o - ^) cos0 O/ 



^jMw-n „ //, an/COS^oA . e r ° ,,Ai/=2\ /cOS0tf 



-f 2 

2— r M 4 Al (fo)(r /o ,V(^))%+ 2^M 2 Al (f )^(x). (86) 

28 



£3 vwv,u,rv~„-,u , £ 3 



This ends the proof of Proposition 14.51 

□ 

Lemma 4.8. With the above notations, for e > small enough, one has 

P x {Ai,2M)(x)) = ^(M 5 Al (fo) -f 2 M 3 Al (f )M*) (87) 

where M^(-) is the non zero entire function defined as M Al (-) + M Al (-). 

Proof of Lemma \4-8[ We proceed similarly as in the proof of Lemma 14.61 Besides remainder 
terms, one must the principal term given by 

I 1 ' 2 {a,ip){x) = p. v.— / — L! — drjipix). 

Using polar coordinates, one gets 
I Al ' 2 (a,ip)(x) 



^ 2 JR2 \V\ £ 2 JR2 l^' 3 



= / / e - 2 e 2,focos e ^ rfr _ ^_ fQ p y / / e -^ e 2,f co S 8 cos() dQ 

^ £ Jr=0 JO E Jr=0 JO r ' 

-2 

where M 3 1 (fo) and M 5 x (fo) are given respectively by (11821) and (1185p . 

□ 

4.1.3 Asymptotic expansion of Ai for 2 < i < 5 

We establish the following proposition for the asymptotic expansion of A, with z = 2, . . . , 5. 
Proposition 4.9. For i = 2, . . . , 5 and e > small enough, one has 

P x (A l (a,^)(x))=0(^). (88) 

Proof of Proposition \4-9[ We proceed similarly as in the proof of Lemma 14.61 
For A 2 (a,ip)(x), we only need to estimate the following term: 



R M {a,4j){x) := (V^(x) + V T ^(x)) p.v. / < ^f L 
By Lemma [4. 4} one gets 



rjdr). 

R2 i??r 



i? A2 (a, V)(x) = ^M 3 Aj (f )(V^(x) + V T ^(x))r/o = O(-). (89) 
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For As(a,ib)(x), we first note that Va^)^ = 0, and 

(n x ,ip(y)) = (n x ,^(x + r]-n x (r])n x )) = (n x ,ip(x) + V^(x)?? + 0(M 2 )) 
= (VjP(x) T n x ,rj) + 0(\rj\ 2 ). 

Thus, we need to estimate the following integral, 

R A3 (a,ip)(x) := % f a em {rj) ^ V ^^ J 1 *' ^ (77 - rjo, rj)dr) n x . 



2 / "s-wv/; I 13 
£ Jr 2 \V\ 



One can clearly apply Lemma [4.31 to R Aa (a, ip)(x) with m — 0,1 and one gets, 



e 2 



and since — = 0(1), one finally deduces that 

€ 



R A *(a^)(x) = 0(±). (90) 



For Ai(a,ilf){x), we only need to estimate the following term: 



R M (a^){x) := ( p.v. / ^^V/, (V^ - V T ^{x))n x )n x . 
Using Lemma [4.41 one gets 



R M M)(x) = ^M 3 Al (r )((V^ - V T ^{x))f)o,n x )n x = O(^). (91) 



For As(a,i/})(x), one gets the estimate 



1 



R A *{a^){x) = 0{-), (92) 



5 



as a consequence of Lemma 14.101 



In summary, for i = 2, . . . , 5, P x (Ai(a, ip)(x)) = 0(— ), which ends the proof of Proposition 



□ 

Lemma 4.10. With the notations above, consider the function defined for x G dQ 



R(x) = r(x,y) -V(aip)(y)da(y), 
Jan 

where r(-, •) is a C|(l) weakly singular kernel and ■ stands for a linear action of r on the 
coefficients 0/ V(cr0)(-). Then, there exists a positive constant Cr such that, for e > small 
enough and x G dQ, one gets 

\\P x (R(x))\\<% (93) 
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Proof of Lemma \4-10[ As done for estimating A\, we use the change of variables introduced 



in Subsection 13.21 and taking into account Lemma 13.41 and Remark 13. 2\ it is easy to see that 
the most 'singular" part corresponds to majorizing 



Va e (r7) 

r B(0,<5) \V\ 

Thus Eq. (193 p follows readily from Lemma 14.31 



drj. 



□ 



4.2 Estimates of the remainder terms W- , i = 2, 3, 4 and j = 1, 2. 

In this subsection, we upper bound the remainder terms P X (W- (x)), i = 2,3,4 and j = 1,2, 
defined respectively in Eqs. ( 158|) . (j5"9l and (J5%|) . More precisely, we prove that 

Proposition 4.11. With the notations above, we have, fore > small enough, j — 1,2 and 
z = 2,3,4, 

P, (W/(*)) = O(i). (94) 

Remark 4.1. One must stress the similarity of our computations with those performed by 

D. Henry in [20] • More precisely, the terms A 1 and A 2 in W^(-), which are (essentially) 

the most 'singular" part in the hypersingular operator E defined in Eq. (1631) . correspond 

to the operator </(•) defined in Theorem 7.4.1 of [20J, page 135, with the specific choice of 
y — x y — x 

Q(x, y, r) = V n (y)- r and n — 3. Also notice that our Lemma 14741 corresponds to an 

\y — x \ \y — x \ 

explicit computation of the polynomial q(-) (cf. Theorem 7.4.1 of [20]) and follows the same 
lines as the strategy proposed in page 137 in [2U]. In particular, one gets from Theorem 7.4.1 
of [20] that W((-) extends uniquely to a continuous operator on dQ. 



Proof of Proposition \4-H\ All the estimates to be established are consequences of (jl73|) - (jl76p 



obtained in Corollary IB. 2 1 We rewrite it as follows. For u of class C 2 and x G dQ, one writes 
47r Eu(x) as the sum of two operators, 



4ttEu(x) = Fu(x) + Lu{x) = p. v. / f(x,y)-Vu(y)da(y)+ / l(x,y) ■ Vu(y)da(y), (95) 

Jdo. Jan 

where " • " stands for an action of the respective kernels which is linear with respect to Vtt(-), 
/(•,•) is a C^(l) kernel (of appropriate matrix size) defined in Appendix IA.2I and the kernel 
/(•, ) defining the singular operator F together with its action is given by 

f(x, y) ■ M(y) := [M(y) + U'iy) ^ J :i + n x n T x [M(y) - M T (y)] * ~_ J , (96) 

for x 7^ y, points on dQ and M is C l matrix-valued function defined on dQ. We are then only 
interested in the first term of the above sum. 

In order to handle the remainder terms P x (W-ys, i = 2, 3, 4, one must handle the evaluation 
at Ki^f of the operators obtained as the composition of Kq defined in (I157p and its iterations 
with W(. In fact, we will show next that all remainder terms P x (W/)'s, % = 2, 3, 4 are 0{\) 
and to proceed, we will be only interested in the contribution of the "most" singular part 
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in each term Wps, i = 2,3,4. For that purpose, we will perform several (and standard) 
reductions. The first one consists in considering the operator Kq instead of Kq since Kq — Kq 
admits a C l kernel. Lemma [4.101 already handles the term W£. Next, recall Kq is a weakly 
singular operator of class C*(l) (see Appendix IA.2I for a definition). To handle the terms W 2 
and Wl, we first need the following result. 

Lemma 4.12. The operator defined on C 1 (9fi) as the composition of Kq and F is a weakly 
singular operator of class C*(l). 

Thanks to the above lemma, the first term in the summation (158]) is controlled as 0(\). 
For the other terms, it is now enough to see that they correspond to compositions of iterates 
of K% with Kq o F and thus we can apply Theorem IA.4I given below on the composition of 
weakly singular operators of class (7) with 7 > 0. We deduce at once that every term 
appearing in the summation (j58|) corresponds to the evaluation at V(V^,-p-)(-) of a weakly 
singular operator of class C% (7), with 7 > 1, and is therefore controlled as 0(^). The term 
W| is handled in a similar way and Proposition 14.111 is established. 

□ 

We now give the proof of Lemma 14.121 



Proof of Lemma \4-lS\ The argument given below is already contained in Section 7.6 of [20] . 



which considers a more general situation (see, more particularly, the proof of Theorem 7.6.3 
page 147, [20]). For sake of clarity, we reproduce the main lines. Let M be a C 1 matrix- valued 
function defined on dQ. Then, the composition (Kq o F)[M](-) is defined, for x € dfl, as the 
sum 

(KqoF)[M](x)=R 1 (x)+R 2 (x), 

where 

Ri(x) = f P-v. / / (X ,~" ,r2( ,f (x-z)(x-zf[M(y) + M T (y)]-^y- da y da z , (97) 
and 



3 ff (x-z,n(z)) 2 (x-z) (n(z),[M(y)-M T (y)}(z-y)) 

R 2(x) = 7- P-v- / / , |4 r— — I , _ |3 da y da z . 

47r JJdnxdQ \z — x\ \z — x\ \z — y\ 

(98) 

Thanks to fll60p . the operator R 2 is clearly more regular than R-y. In the sequel, we only 
provide details for R\ and only give the estimate for R 2 . 

We next develop in coordinates the above expressions and obtain that, for i — 1, 2, 3, 

3 



(R 1 (x)) = j^-X>v. / (M(y)) kl da y 
V ^ 47r ^ Jan 



(x - z,n(z))(x - z)i(x - z) k (z-y)i 



an 1 \ z-x\ 5 \ z - y | 3 



| (x-z, n(2:))(x - - z)i (z - y) k 
\ z — x \ 5 \ z — y \ 3 



da z . (99) 



32 



The integrand of fl99|) shows that i?i is the contraction of M(-) and a tensor field of order 
(2, 1) defined (in coordinates) by the interior integral in fl99|) . In order to describe R\ as a 
convolution, we prefer to rewrite fl99l) in a more elementary way, as follows, 



(Rx{x)) = -f-p.v. / Tr(M(y) 



Ci(x,y)) da 



where the kernel Ci(x, y) is defined for x ^ y and 1 < i < 3, as 

Ci {x,y) :=p.v. / 1 : I is rr— - is — + 1 y _ Tf 13 — d(Jz - ( 10 °) 



— 2/1 I z ~ y 



Let (ej)i<j<3 be the canonical basis of 1R 3 . Then one has Ci(x,y) = di(x,y) + di(x,y) T , where 

di(x,y):=-p.v. k (x,z)[g\z,y)\ da Zl (101) 
Jan 

i.e., di(x,y) is the kernel corresponding to the convolution of Kq with kernel k°(-, •) given by 
k°(x,y) :-- 



1 (x-y,n(y)) (x - y) (x - y) T 



x — y | | X — ?/ | " | a; — jy | | x — y 
and the singular operator G l with kernel (/(•,•) given by 

e;.(x - y) T 



9 l {x,y) :-- 



x-y 



To perform that analysis, one writes fllOip in the chart h x defined in (|46j) and only considers 
the most "singular" term of the composition, which is given by 

p.V. f ^ m Ti^l 2) 

Jb(o,s) \vr Iv-Vvr 

Here, ei is the orthogonal projection of onto T x dQ. In f ll02p . one clearly recognizes the 

convolution between the kernels — — r^VV T an d % r? - The first kernel can also be written as 

I V I 1^1 

l^|Q(|^|) where the components of Q are homogeneous polynomials of degree four defined on 
S 1 . According to [201 Th. 7.3.1 p. 128] (which refers to [32] for more complete computations), 
the Fourier transforms of these kernels are respectively equal to 



and 



where 71 is a positive constant and the components of Q are homogeneous polynomials of 
degree four. We get that the Fourier transform of the operator whose kernel is given by (11021) 
is equal the product of the two Fourier transforms written previously and, as a consequence, 
that operator is weakly singular of class C 3 (l). The same conclusion holds true as well for R\. 
A similar line of reasoning shows that R2 is weakly singular of class C 3 (2) and Lemma [4.101 is 
finally proved. 

□ 
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5 Proof of Theorem 



1.1 



In this section, we establish in full generality the Foias-Saut conjecture in 3D as stated in 
[IH]. First of all, notice that there is a countable number of resonance relations as defined in 
Definition 11.11 To see that, simply remark that, for every positive integer N, there exists a 
finite number of resonance relations of the type A& = Y^j=i m o^h "with Ai < ■ ■ ■ < A^ < A^, so 
that k + Y^j=i m j — We use (RR) n , n > 1, to denote these resonances relations. 
Fix a domain Qq £ with £ > 5. We define, for n G N, the sets 

A := O?(fi ), 

and, for n > 1, 

A n := {Q + u, ue W^°°(fi , M 3 ), Q + u G A 
and the n first resonance relations (RR)j, 1 < j < n, are not satisfied}. 

Set A := f] l&N A n . Note that 

A = {Q + u, u G W^ £ '°°(fio,ffi 3 ), fi + « G A (SX>n 0+u ) is not resonant}. 

For n > 0, each set is open and one must show that A n+ i is dense in A n . Reasoning by 
contradiction, assume that there exists n G N so that A n+ i is not dense in A n and fix (RR) n+ i 
to be equal to A& = Y^j=i m j^ji f° r some integers k, I, mi ■ ■ ■ , m;. With no loss of generality, 
we assume that there exists fl G Of and e > so that, for u G W e '°° with IHI^oo < £ 5 we have 

(i) the /c first eigenvalues X\(u), . . . , Afe(n) of (S-Dn+u) are simple; 

(ii) the resonance condition holds true: 

i 

X k (u) = J2 m Mu)- ( 103 ) 



By Condition (i) and Eq. (I3"3"|) . one has, for u G with ||if||£ i0 o < £ and 1 < j < k, 



A» = - /J^ll^lh (104) 



where 0j is the orthonormal eigenfunction associated to the eigenvalue Xj of (SDq). 
Taking the shape derivative of Eq. (11031) , we have 



Since Eq. (I105[) holds true for all u small enough, we obtain 



^|| 2 -V mj |ft|f = fl on dQ. (106) 
On ^ On 
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Continuing the argument by contradiction, we assume that Eq. (I106p holds true in a 
neighborhood of Q and take again the shape derivative. By Proposition 13. 2[ we have, on dQ, 



dv ' dn 



-E 

3=1 



m 



dv ' dn 



(u,n) ( 



d d(f) k 



dn ON ' dn 



>-E 



5 



' <9n <9iV ' dn 



) • (107) 



We choose a variation it such that (it, n) = V n with V n defined in Section 13 .21 Using 
Proposition 13.51 together with Eq. ( llOSp . since fj is an arbitrary unitary vector of IR 2 , we 
obtain, on dQ, 



T 



dn V dn 

From now on, fix x G dQ such that 



E 



<9n V dn 



0. 



108) 



<9n 



(x) 7^ 0. Recall that such an x exists by the result 



of Osses in [27]. According to Eq. fll08p . there exists an open neighborhood O x of x on dQ 
such that, for 1 < j < I, there is a C 2 function fij such that 



dn 



dn 



on O x . 



(109) 



In addition, one has, 

i 

1 - m jV 2 j = 0> on o x . (110) 

It is clear that all the equations from f)104p to (IllOp were obtain by assuming that Eq. f ll03p 
holds true in an open neighborhood of u = 0. As a consequence, these equations must also 
hold true in an open neighborhood of u = and thus, one can take the shape derivatives of 
Equations fllQ8p at u — along any variation. We will perform such a shape derivation along 
the variations V n defined in Section 13.2} with this time the real number 5 > chosen so that 
the support of V n is contained in O x . Using Lemma I3~TT| the shape derivative of Eq. f)108p is 
equal to 



dv V dn 



dn V dv 



3=1 



.1 



Pk 



dn 



n) 



n 



dn 



T 



dn 



k T 

n 



dv V dn 
i 



+ 



E 

3=1 



dn V dv 



,n) 



T- 



dn 



n 



dn 



d 



dn dN V dn 



d(j) k r d d(j) k \ T _ 
dn \dndNJ ^ 

3=1 



d 



dn dN V dn 



'111] 



dn 



d 



3 T 

n 



•V 

dn V dn dN / . 



where the above equation holds on dQ. 
Moreover, on O r , one deduces that 



d d(j)j 
dn dn 



V(/i j (V0 fc n))n 



dnj 



dn dn 



+ fijV 2 (f)k(n,n). 



h V T (f) k n + /i j V(V T fc n)n = /i j V(V T 0A,.n)n. 



V(/Xj(V%n))n 



dn 
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This implies that, on O a 







dfij 



+ (J>jVk, 



dn dN dn dn 
with Vk '■= V 2 0fc(n, n) + V(V T 0fcn)n. Therefore, one has on O x 



(112) 



d d(j) k 
dn dN \ dn 

i 



+ 



(j^(tyk\ T r d d<t>j / 

<9n V «9n <9iV / ^ m ^ j Idn dN \ 
r 



i=i 



<9n 



<9n 



A; T 



3=1 



<9n diV V dn 
dfi 



+ 



9 d(j)j\ 



T- 



dn \ dn dN ) . 



dn dn V dn 



3=1 



dn dn V dn 



Plugging Eqs. ffT09|) . flTTU]) . and ( TTT3]) into Eq. ffTTT]) . we obtain on O s that 



113) 



<9i^ V <9n 



+ 



<9n V dv 



3=1 



dv V dn 



T 



+ 



<9n V dv 



i 



Vk - Yl "'ji'ji'j 

3=1 



n 



dn 



dn 



k T 

n 



i 



2V n Y m j^j 

3=1 



dfij 



dn dn V dn 



On O x , set 



3=1 



Wk_ d Ay 

3 dv dv 



(114) 



(115) 



The main part of the rest of the argument consists in deriving the main term of the 

asymptotic expansion of P x (dk) at x, in terms of the powers of -. The first step will be 

e 

to establish that P x (dk(x)) = O(-) and, in the second step, we will compute precisely the 
coefficient d l defined as 

P x (d k (x)) = j + 0(1), 

where the coefficient d 1 will depend on the parameters involved in the special variation V n . 
Once this is performed, we will resume the contradiction argument using the information 
contained in d 1 . 

To prepare these computations, we first rewrite Eq. (I114p using again Eq. (IllOp as 



3=1 



dn dn \ dn 



3=1 

T 



dn 



dn 



k T 

n 



(116) 
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Multiplying Eq. (11161) from the left by ( — — ) and from the right by — — , we obtain the 

v On ' On 
following scalar equation which will be used to achieve a contradiction. 

(^, 4 >=[£m,„,^]K. (117) 

We now prove the following lemma. 
Lemma 5.1. With the notations above, one has P x {dk{x)) = 0{—). 

Proof of Lemma \5. 1\ Set ip := and, for y £ O x , 

on 

i 

P(v) := ^2 m iN( x )(N(v) - H( x ))- ( 118 ) 

Then, one has (3(y) = 0(\x — y\ 2 )- More precisely, if we use the parameterization defined in 
Eq. (14"6)) . we obtain (in local coordinates) 

0{V) = \ (^(x) V T K xV + r, T H x ri) + 0(M 3 ) := \n T F x n + 0(\ V \ 3 ), (119) 

where H x denotes the Hessian matrix of (3 at x. Note that by taking twice tangent derivatives 
of Eq. flllOp . we know that H x is a negative semi-definite matrix. 

Consider now the representation formula of as described in Eq. (j!62p . Note that, in 

P x (dk), two contributions give rise to the term of order of O(-), namely these coming from 

and respectively. 

The term corresponding to in that equation is equal to 

P x (E(fivjy±)(x)). (120) 

Thanks to the estimate of (3 in (11191) . it is clear, by proceeding as in Subsection I4.2[ 
that all the other terms of the representation formula of P x (dk) are indeed of the type O(-). 

Therefore, one has only to determine the asymptotic expansion of the term given in Eq. (I120|) . 
According to Lemma [4.2[ it amounts now to estimate the five terms P x (Ai), 1 < % < 5, and 
after elementary or standard computations using systematically Eq. f 1 1 1 9 j) , we obtain 

P x {Ai{a,W){x)) = ^^ 2 ^r / T F^((^(x),r / )(r / -r /0 ) + (r/-r /0 ,r / )^(x)) C /r / +O(l), (121) 

P x (A 2 {a,Pif>){x)) = ~ [ ^U(x)n T F x n+^(x),n)F x n)dn + 0(l), (122) 

and, for 3 < j < 5, A 5 (a, Pifj) {x) = 0(1). 

Let us now treat the terms given bye (A) . Note that the presence of these terms reflects the 
fact that (pj and 4>k correspond to different eigenvalues of the Stokes operator. Using Lemma 
IA.1[ we define the operator A x (a,ip) as follows 

1 A ' 8nJ R2 \rj\ \rj\ 1 K J 
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It is clear that 

A A («,^)(x) = 0(i). 
Therefore, with above notations, we have 

i 

P x (d k (x)) = A x (a, ip)(x) + A 2 (a, i/j)(x) + ^ m j fi j (fi j A Xh (a, i/j)(x) - A x > (a, fijip)(x)) + 0(1) 

3=1 

= 0(- £ ). (124) 
This ends the proof of Lemma 15. 1[ 

□ 

Let us now pursue the proof of Theorem ll.lt Since the value of the right-hand side of Eq. 
(I117p at x is given by the following expression 

/ -2 

E^,(x)^(x)]^ = 0(l), 



we conclude that 



which implies that 



£m^-(s)^(s)=0, i.e., f(*) = 0, (125) 



, -(x),d k (x)) = 0, (126) 



dn 



In order to get additional information from Eq. ( I126p . we compute explicitly the numerical 
coefficient in front of - in the asymptotic expansion of P x (d k (x)). It is enough to have a closer 

look at the representation formula of P x (dk) as described in Eq. ( I162p . From Eq. (jl24j) . we 
have 

P x (d k (x)) = a 1 + a 2 + pa 3 + 0(l), (127) 



where 



f ( 1 

&1 := 4vrc2 J 2 'J^^FxvfaWirfiv-Vo) + (rj-Tjo,Tj)iJ}(x)Jdi], (128) 

a 2 := —1- f ^j^U{x) V T F xV +(^{x), V )F xV )d V , (129) 

1 f a £ (r]){jj(x),r]/\r]\) r) 

a 3 := -— / j-j — -dr), (130) 

8tt J R 2 \r]\ \rj\ 



i 



p ■= J2mj^(X k - Xj). (131) 

i=i 

Notice that p > since > Xj for 1 < j < I and at least one of the integers rrij is positive. 
We now compute the coefficients a iy 1 < i < 3. For 9o £ S 1 , we set 

_ fcos9 -sin# \ ._ R Tpn — fp*^ 

^ - ^sin^ cos^ J ' Fe ° - ^o^^o - (Fe hj=i,2- 
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For 1 < i < 10, we define the functions Mi as follows. 



M 6 {z) := / e~ r \ 2 dr / e 2rzcos9 d9, (132) 
Jo Jo 

pOO t>2n 

M 7 (z) := / e~ r \ 2 dr / e 2r z cos 9 cos 2 9 d9, (133) 
Jo Jo 

poo 

M 8 (z) := / e~ T \ 2 dr \ e 2r z cos 9 cos 4 9 d9, (134) 
Jo Jo 

poo p2n 

M 9 (z) := / e~ T \dr \ e 2r z cos 9 cos 9 d9, (135) 
Jo Jo 

poo p2n 

M 10 (z) := / e~ r \dr \ e 2rzcos9 cos 3 9d9. (136) 
Jo Jo 

The expressions of a\, a2, and 03 are summarized in the following lemma whose proof is 
postponed in Section [RBI of Appendix. 

Lemma 5.2. We have 



1 e" f o 



a 1 



Air e 

[2F 22 M 6 (f ) + (2F£ - 3F 22 )M 7 (f ) - (F# - F 22 )(M 8 (f ) + M 10 (f )) - F 22 M 9 {f )]^{x) 
- m 2 M 9 (f ) + (FH - F 22 )(M 10 (f ) - 2M 8 (f )) + (^ - 3F 22 )M 7 (f ) 
+F 22 M 6 (f ))(^(x),r /0 ) + 2F e 1 o 2 (M 9 (f ) - M 10 (f ))<^), 

2F, 12 (M 7 (f ) - M 8 (f ))^(x) x +4F e 1 o 2 (M 7 (f ) - M 8 (f ))^(z),%)^), (137) 



_> 



« 2 = -^^{[i^M^C^ + CFj-Pg'^Cfo)]^) 

+F E ((M 5 Al (ro) - M 1 Al (f ))^(a:) + (2M^(f ) - M 5 Al (f ))(^(x), fjo)fh) }, (138) 

(M 5 Al (f ) - M^(f ))^(x) + (2M^(f ) - M 5 Al (ro))(^(^)^o)d • (139) 



1 e" f o 



«3 



in e 



Let us now finish the proof of Theorem 11.11 We choose 7/0 L ^(x). Without loss of 
generality, we also assume that fj = (1 ; 0) T and ^(x)/|^(x)| = (0, 1) T . Recall that we have 
chosen x such that ip(x) 7^ 0. Then, we deduce from Eq. (I127p and Lemma [5.21 that 



d k {x) = ^ (axip(x) + a 2 ip(x) ± + a 3 F x ^(x)) , (140) 

where, 

ai := 2F 22 M 6 (f ) + (2F x n - 3F x 22 )M 7 (f ) - (F, 11 - F 22 )(M 8 (f ) + M 10 (f )) - F 22 M 9 (f ) 

-((Ff + \ )M*(f ) + (Fi 1 - F 22 - |)M*(f )), 
a 2 := 2F a ! 2 (M 9 (f ) + M 8 (fo)-M 10 (fo)-M 7 (f )), 
« 3 := -(M^^-Mf^fo)). 
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Plugging Eq. (I140p into Eq. (jl26j) . we obtain 

ai + Ff a 3 = 0. (141) 

The final contradiction will be obtained by showing that the two non zero entire functions 
of fo given by ot\ and «2 cannot satisfy Eq. ( (11411) ). For that purpose, we need to get more 
explicit expressions. Eq. ( (I14ip ) writes 

= 2F? M 6 (fo) + {2F1 1 - 3F? )M r (f ) - (F x n - F? )(M 8 (f ) + M 10 (f )) - Ff M 9 (f ) 
-((2Ff + ^)M^(f ) + (F^ 1 - 2Ff - |)M*(f )) 
= F x n (2M 7 (f ) - M 8 (f ) - Mio(fo) - M^(f )) 

+ Ff (2M 6 (f ) - 3M 7 (f ) + M 8 (f ) + M 10 (f ) - M 9 (f ) - 2M 5 Al (ro) + 2M 1 Al (f )) 

+ £(M*(f )-M*(f )). (142) 

Since the right-hand side of Eq. (11421) is an entire function, we deduce that all the coef- 
ficients in its power series expansion are equal to zero. We need the following lemma whose 
proof is deferred in Section IB. Tl of Appendix. 

Lemma 5.3. The entire functions involved in Eq. (11421) have the following power series 
expansions 

2M 7 (z) - M s (z) - M 10 (z) - M^(z) 
= E (2^)! ^ + 2 )/2 ^ 2 2p + 3 Z - g (27TT)! + 2 } 2 ' 



2M 6 (z) - 3M 7 (z) + M 8 (z) + M 10 (z) - M 9 (z) - 2M^(z) + 2Mf 1 (z) 
^2^+1, l w 6p 2 + 16p+| 9n ~ 2 2 ?+ 2 w 3 W 1 



o J2pjl iP + 2 ) ^(2p + 2)(2 P + 4) 2 - E (27TI)! + 2 )/2p+2 2p + 4 



2p+l 



Using Lemma I5T31 and considering the coefficients of the odd powers of z in the power series 
expansion of the right-hand side of Eq. (I142j) . we deduce that 

Fl l (2p + 3) + Ff = 0, for all pel. (143) 

This implies that 

Fl l = Ff = 0. (144) 
Therefore, using Eq. (I142p . we have 

£(M*(f o )-M*(f o ))=0. (145) 

Recall that M^ 1 (z) — M^ 1 (z) is equal to M^ 1 (z), then not identically equal to zero. Then 
Eq. ffl45D yields that 

p = 0, 

which is in contradiction with the fact that 

i 

P = E m il j2 Mk - Aj) > 0. 
J'=l 

Theorem 11.11 is finally proved. 
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A Layer potentials and representation formulas 



r 



Most of the material presented here is borrowed from [H [3j [231 120] • Let A be a non negative 
real number. Consider </> and p satisfying the eigenvalue problem associated to the following 
Stokes system 

(A + A)</>-Vp = h mQ 
div (j) = in Q 
<p = g on dfl 

V = 0, 

n 

under the compatibility condition 

■ n ds = 0, (146) 

where n is the outward unit normal to dQ. Recall that, for such a pair of fields, the conormal 

derivative denoted by — — was defined in (fTTI) . 

ov 

A.l Layer potentials 

We denote by <9j the operator and by \J — A the complex number iy/X. In this section, 
we adopt the Einstein summation convention which omits the summation sign for the indices 
appearing twice. 

Fundamental tensors We define the fundamental tensors T x = (r^-)| - =1 and F = (-Fi) 3 =1 
as 

x.. P \T->\x\ i pV-^kl _ i 



h3 4ir\x\ 4vrA 4 3 



(147) 



FAX 



4ir\x\ 3 

In the sense of distributions, straightforward computations of the fundamental solution of 
Helmholtz operator A + A allow to get 

(A + A)r* - djFi = S^x), and = 0, 

where we use S(x) to denote the delta distribution based at x G M. 3 . The tensor T , which is 
the fundamental tensor for the standard Stokes system, is defined as 

Ti(x) : = 



1. 


( % _ 


- 


8tt 


V \ x \ 


\ x 3 



and one has, uniformly on compact subsets of 1R 3 , 



r*(x) = i%(x)-*^ + o(\). (148) 



Denoting that 
we have 



A* (*):=!* (*)-!*(*), (149) 



Sij v — A A 



A£(*) = " 32^ A ^(^) + 0(|x| a ). (150) 
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where Ay(-) is defined by 

A tJ (x):=35 i:i \x\-^-. (151) 

After simple computations, one gets the following useful result. 
Lemma A.l. We have 



d 2 A x (x-y) A i(n x ,n y ) (x - y)(x - y) T n y n T x 



dN(x)dN(y) 8vr 
where T\ is a kernel of class C 1 . 



+ 



\x — y\ \x — y\ 2 \x — y\ 



+ T\, (152) 



Single and double boundary layers In the sequel, we use the Einstein convention 
for summation signs, i.e., we omit them for indices appearing twice. Let <fi = ((ft 1 , <fr 2 , <f) 3 ) e 
L 2 (dQ) 3 . The single-layer potential pair (5q, J-q) with density <fi is defined, for x G Q, as 



S^Ux) = / r*(x-y)(P(y)da y , 1 < z < 3, 
'an 



(153) 



^nM(a;) = / Fj(x-y)(j) J (y) da y , 

IdQ 

A 



while the double hydrodynamic potential pair (Dq,Vq) with density is defined by 



Ian ( < 



D&Mz) = I \^^(x-y) + F l (x-y)n 3 (y)j ( p j (y)da y , 1 < i < 3, ^ 
/" <9F 

V n [0](x) = -2/ —L(x-yW(y)n l (y)da y . 
Jdn ox i 



We quote from [23] that 



dN(y) y dyi dyj 

Some background results about the layer potential representations From [I], we 
quote the following integral equations satisfied by A and the associated pressure p x . First, 
we have the following representation formulas, 

<P\x) = -S^}(x) + D^](x), xGfl, 

j&\ ( 155 ) 
p\x) = -^[^-](x) + V n [0 A ](x), X6S1. 

Applying the trace stress operators and taking into account the single layer potential as well as 
the jump relations for the double layer potential across the boundary, we get for belonging 
to L 2 (dQ) 3 the following relations, 

DM(x) = (h + K^Mx), a.e.onffi, 
a 2 X (156) 

— SM(x) = (--I + (K>y)[cf>](x), a.e.on9fi, 
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where the kernel -KqM * s defined a.e. on dQ by its components, 



K&tWx) :=p.v. 



dT ij 



(x-y)(j>>(y) da y + p.v. / F^x - y)<jp (y)nj(y) da y . (157) 



Here, the notation " p. v." indicates the Cauchy principal value when the integrand is singular 
at x, more precisely 



p.v. 



lim 

e-s-0 



>n\B(x,e) 

where B(x, e) is the ball centered at x of radius e. The adjoint operator Kq* of Kq is defined 
similarly a.e. on dfl by its components 

r dr x 



^*Mi(x) = p.v. 



(x - yW{y) da y - p.v. / Fi(x - y)(f) J (y)n j (x) da y , (158) 



for all functions belonging to L 2 (<9f2) 3 . Let us recall that in the case of the standard Stokes 
system (A = 0), we have 



47T 



I" / (a? -2/) 



[x - y,n(y)) (x - y,(f>(y)) 
\ x — y \ 5 



d<Jy. 



(159) 



An important fact is that the single and double layer potentials Sq and are compact 
perturbations of the single and double layer potentials corresponding to the standard Stokes 
problem. 

From the C regularity of the boundary V with i > 4, it comes that 



\(x - y,<f)(y))\ < C\x - y\ 



(160) 



hence, we deduce (cf. [23]) that the mapping K^cf)] : C a (dQ) (-)■ C a+1 (dQ) is in fact con- 
tinuous. That shows that -K^[</>] has a weakly singular kernel and then that it is a compact 
operator on L 2 (dQ) 3 . According to (I148p . the operators — and — Dq are smoothing 
operators. 

Thanks to the integral representations provided in the preceding paragraph, we can use 
the trace and the stress operators to deduce the second boundary integral equation satisfied 
by the conormal derivative. Indeed, by using the same arguments of jump relations and the 
integral equations satisfied by A , we get 



dv 



\x 



9DM 



dv 



x 



p.v. 
p.v. 



d 2 T^(x-y) 

m dN(x)dN(y) 
d 2 T%{x-y) 

m dN(x)dN(y) 



<fi J (y) da y 
(jy'iy) da y 



(161) 



+ 



d 2 A*(x-y) 
m dN(x)dN(y) 



^(y) da y . 



We cannot deduce directly the Neumann data (conormal derivative) since the operator y—1 + 

(Kq)*j is not invertible. We give, in the next paragraph, the recipes to get the solution of 
the system by using the projector methods. 
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A. 2 Weakly singular integral operators of exponent a > 



The rest of the paragraph follows Section 7.2 of [20]. Recall that the conormal derivative is 
solution of Tx = f where T = I + 2(Kq)* is a Fredholm operator with a nontrivial kernel. 
We use 1Z(T) to denote its closed image and Af(T) its finite dimensional null space. We can 
therefore find projections P and Q of finite rank such that there exists a unique operator S 
satisfying TS = I — Q and PS = 0. Hence, the equation Tx = f has a solution if and only if 
Qf = 0. In our context, we have T = I — C with C = —2(Kq)*, which is a compact operator. 
To proceed, we need some regularity assumptions on the operator T. For that purpose, we 
recall the following definition [201 Definition 7.1.1, pi 17]. 

Definition A.l. Let A be an open set in M 3 . A function K(x,y) defined for x ^ y in A x A 
is a kernel of class Cl(a) in A (r non negative integer, and a > 0) if it is C r for x ^ y and for 
any S > and \i\ + \j\ + \k\ < r, one has 

di& y (d x + d y ) k K(x, y) = 0(1 + \x- y \<*-™-W-\i\-6)^ 

uniformly for x ^ y in compact subsets of A. If a > m + \i\ + we require d % x d 3 y (d x + 
d y ) k K(x,y) to extend continuously to {x = y}. 

Assume now that T = I — C is an integral operator where C has a kernel belonging to 
Cl(a) for some a > 0. We may choose the projections P and Q to be integral operators with 
C r kernels so that, if S is the operator such that 

TS = I-Q, 
S = I + R, 
PS = 0, 

then the resolvent kernel R is an integral operator with CJ(a) kernel. In fact, R — (C + C 2 + 
• • • + C J ) has kernel of class C£((j + l)a) for each j > 1. Hence, for A r sufficiently large, the 
operator R — Ylk=i ^ nas a smo °th kernel of class C r [281 Chapter III, p 79-94]. In summary, 
one has the following result. 

Theorem A. 2 (Theorem 7.2.3, page 125 in [20]). We suppose Q regular of class C r+l , for some 
r > 0. If K a kernel of class Cl(a), then we may choose the kernels P and Q of the projections 
to be of class C r and such that the resolvent kernel R belongs to Cl(a). Furthermore, the kernel 
of R — (K + K 2 + . . . K ) is of class C r for N large enough. 

We return to the study of Eq. (116ip . We introduce the vectors b^\x) = (6^ ' l (x))i and 
e^ x \x) = (e[ X \x)) where 



and where 



Hence it comes that 



(A) i , / - y) , 



N N N 



6 (o) + (^ir fe )6(°) + (^K fe )eW + (i2 - Y,K k )(b^ + e W). (162) 

fe=l fc=0 fc=l 

With our specific choice of Dirichlet data, we will show that N — 1 is sufficient in our 



<9^ . 
spec 

context and that all the other terms in the sum will be absorbed by the remainder. 
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A. 3 Composition of weakly singular kernels 



For applications to our result on generic perturbation of the boundary, we need to give an 
explicit representation of the conormal derivative or at least, of its principal and subprincipal 
parts as it is treated in the case of the Laplacian (for more details in the Laplacian case, one 
can refer to [31] ) . Some preliminaries are required in order to study the resolvent kernels and 
their regularity. We begin by recalling some results due to D. Henry (cf. [20] ). It concerns 
kernels K(x, y) of the form 



where Q(x, y, s) is of class C r (r > 0) on IR 2 . We will denote by /C(r, a) the set of such kernels, 
which is a subclass of CJ(a). 

These kernels are in fact smoothing operators and we recall the main result of [20] . 

Theorem A. 3 (Theorem 7.1.2 in [20J). Given a kernel K belonging to the class JC(a,r), 
a, r > 0, we denote by K the corresponding integral operator 



Then we have 

• K : W j ' p !->■ W k,p is a compact operator if j — ja > k — —; 

• K : C^ a C k,T is a compact operator ifj + a + a>k + r, k < r and k < j + a. 

As it was mentioned in [20], the above result can be summarized by the fact the operator 
K is smoothing of order a. We will also need a result on the composition of certain weakly 
singular operators. For that purpose, we first define the composition of corresponding kernels 
as follows. 

Definition A. 2. Let K and L be kernels belonging to K,(a,r) and K,(/3,r) respectively with 
a, /3, r > 0. Then K o L is defined by 



Then, one has the following property. 

Theorem A. 4 (Theorem 7.1.3, p. 119 in [20J). Let K and L be kernels belonging to )C(a,r) 
and /C(/3,r) respectively, with a,(3,r > 0. Then K o L is kernel of compact support belonging 
to K(ct + (3, r). Furthermore, ifa + (3>r + 2, then K o L is of class C r . 



To these kernels, are associated integral operators u h-> / K(x, y) dS(y) where dS is the 



surface area measure on dQ. In a first step, we begin to work in IR 2 . To transfer all the results 
to dfl (in particular, those provided above), one has to follow the classical steps: construct a 
partition of unity and then define the integral by a local change of variables as it is precisely 
performed in [20J, Section 7.1]. 




(163) 





(164) 
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B Proofs of computational lemmas 



B.l Proof of Lemma [278 



From Eq. (12811 . we get the following system 

- (A + A)#(u) + Vp'i(u) = Ki u )Hu) infi, (165) 
div #(u) = in n, (166) 

^(u) + (u-n)^^ = ondQ, (167) 
an 

p'iiu) + div (upi(u)) G Ljj(fi). (168) 

Multiplying ( 116511 by <p k (u) with 1 < k < m, integrating over Q and using Corollary 12. 3\ we 
have 

Al(«)fe = - / M«)[(A + A)#(«)-VpJ(«)]= / #(«)-^ 



Hence, it comes that 

\m k = ~ jju-n)—.^-. (169) 



Moreover, by Lemma (12. 4p . we have 

d(j>i(u) d(j) k {u) d(pi(u) 



dn dv dn dn 



+ V T 4>k{u)n- p k (u)n) 



u) d(j) k (u) d<t>i{u) ,dcj) k {u) T T dcj)i{u) d(p k {u) 

n ) n - 



dn dn dn dn dn dn 

Therefore, we immediately get Eq. 



B.2 Proof of Lemma 14.2 

Lemma [4.21 is derived from [2"Tj Lemma 2.2.3 Formula (2.2.34) and Lemma 2.3.1] by straight- 
forward computations. For the reader's convenience, we first summarize these results in the 
following lemma and then give the proof of Lemma 14.21 

Lemma B.l. Let dQ be of class C 1 andu = (w )^=i,2,3 be a Holder continuously differentiable 
function. Then the operator E defined in (l62|) can be expressed as follows 

1 f 1 

Eu(x) = - — (n x x V*) • / , An y x V y )u{y)da{y) (170) 

47r Jan F — U\ 

-^M(d x ,n x ) [ ix ~ y){x ~ y)T M(d y ,ny)u(y)da(y)) (171) 
An Jan F ~ V\ 

+-r(y2rnik(d x ,n x ) r(m kj (d y ,n y )u e )(y)da(y)) , (172) 

47r \ fc=1 Jan F ~~ y\ /j=i,2,3 

where the £ th -column of the matrix (n y x V y )u(y) is given by the vector n y x V y u (y), and the 
Giinter derivatives M. is given by the following matrix of differential operators 

M(d x ,n x ) = (m jk {d x ,n x ))j, k =i,2,3 ■= (n x , k d Xj - n xd d Xk ) jtk =i,2,3, 

with n x = (n XJ ) i= i j2 ,3- 
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Corollary B.2. Under the assumptions of Lemma \B.l\ we have 

4ttEu(x) = p.v. / p^Ahu(y) + V T u(y))(x-y)da y (173) 
Jan F ~~ V\ v y 

+ P.v.y^ ^ (174) 

- £ V- ^ ( Vn(y) + vTu ^ y) ) nJa y (175) 

/" (x-y,n y )f Q (x-y)(x-y) T \ 

Proof of Corollary \B.t& For (I170j) . we get 

(n* x Va.) • / , \ , (% x V y u\y))da{y) 



= p.v. / (n^V^ — ,) • (rij, x V/(t/))(f(x(?/). 
Jan F — 3/1 

For x ^ y, one has 



(n* x V x - ) ■ {n y x V y u\y)) 
\x y\ 

{n T x n y ){V x -^—V T y u\y)) - {V y u\y)n x ){V x —— 
\x — y\ \x — y\ 

n T x n y . (x-y) T n 



hi 



V y u\y){x -y) + — ^V/^n 



\x — y\ 3 \x — y\ 



Therefore, we have 



[n x x V x ) • / | — - — An y x V y )u(y)da(y) 
Jon F — V\ 

-p.v. / p^V y u(y)(x-y)da y + p.v. [ {x . - V ' ^ V y u{y)n x da y . (177) 
Jan F ~~ £/l -/an F ~~ 1/1 
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We compute now the second piece of (11711) and obtain for x ^ y 

M{n x ,d x ) : - = y m ik (n y ,dy) , J - J — 

\x - y\ s \^ \x — y\ 6 /i,i=i,2,3 



(x k - Vk)(xj - yj) 



Xk \x~y\ 3 A,j=l,2,3 

3 



E n r 3 (xi-yi)( x k-yk)(xj -yj) + Sikjxj - yj) + 5ij(x k - y k ) ] 



k=l 



x — y\ 3 \x — y\ 3 \x — y\- 



_ n r p ( x k - yk) 2 {xj - yj) | [xj - yj) 5 kj {x k - y k ) . 

\x-y\ 3 \x-y\ 3 \x-y\ 3 s Jij=i,2,3 

z (n x ,x-y) _ y) ,_ y) T + n x (x-y) T + (n x ,x-y) ^ 

\x — y\ 5 \x — y\ 3 \x — y\ 3 

3 n x (x-y) T _ ^ n x (x - y) T _ n x (x - y) T 
\x — y\ 3 \x — y\ 3 \x — y\ 3 

(n x ,x-y)( Q (x-y)(x-y) T 



\x — y\ 3 \ \x — y\ 2 

Therefore, we have 

M(d x ,n x ) [ (X ~ V){X ~ V)T M(d y , n y )u{y)da{y) 
Jan \x-y\ 3 

f (n x ,x-y) ( [x -y)[x -y) T \ . , . , , . 

= P,V - J m \ x - y\ 3 \ 3 ~ \x~^yj 2 )M(d y ,n y )u(y)da y , (178) 

keeping in mind that there is no principal value if one uses (11601) . 
We finally turn to (11721) . One has, for x ^ y, 

3 - 

Y~] (mik(d x , n x )- r)(m kj (d y , n y )u\y))) 

rr -1 , \ x ~ y\ /j=i>2,3 



!,fc=l 
3 



.'/I 

i,k=\ \ y\ \ y\ > ' 

("x,(V«(?/)-V T tt(t/))(a;-?/)) (x-y,n y ) T 

~ \*-y\ 3 u{v)nx 

+ P^V T u(y)(x-y). 



\x - y\ 3 
Therefore, we have 

3 



( y2 rni k (d x ,n x ) [ - — - — :(m kj (d y ,n y )u e )(y)da(y)) 

V f k ~Z 1 Jdn\ x -y\ '.7=1,2,3 



n x ,(Vu(y)- V T u{y)){x-y)) f (x-y,n y ) T 
■ riyda - p.v. / — -f-V u(y)n x da y 

\ x ~y\ 3 Jan \ x ~y\ 3 

(n x , n y ) nT _ 



+ p.v. / ^^LLv 1 u(y)(x-y)da y . (179) 
Jan \ x ~ V\ 
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Gathering (1177)1 . ffTTHj) . and (fTT9]) . Corollary [Bj is proved. 

□ 

Proof of Lemma \4-2\ Recall that u = aip with a : dQ h-> R and ^ : dfl i— >• R 3 . We note that 

V(aV>) = aV^ + ^Va, . (180) 

and 



1,2,3 



2.3 



M(d y ,n y )(aip)(y) = (^2m ik (n y ,d y )(a(y)ip k (y)))_ 

k=l % ~ 
3 

= a(y)M{d y ,n y )i)(y) + ( ^K,^a(i/) - n y>i d yk a(y))il) k {y) J 

fe=i 1 
= a(2/)A1(9 y , ny)if>(y) + ^(y))V T a(y) - {Va{y)^{y))n y 
= a{y)M{d y ,n y )xl){y)-{Va{y)i){y))n y . (181) 

Then, the expressions of Aj, 1 < i < 4, simply result from developping Vw in (11 73)) and (11741) 
of Corollary IB .21 and A 5 collects (1175)) and (11761) as a weakly singular operator of class C 3 (l). 
Hence Lemma [4.21 follows. 

□ 

B.3 Proof of Lemma 14.31 

Using polar coordinates, we have 

JB(0,<5) I ?7 I £ </b(0,5) I 7 /! ^ Jo JO £ £ 



e o 



£ l-m 



<5/e /-27T 

exp (— r 2 + 2rf cos 6)r rn drd6 



o </o 

2 "OO /"27T 



< -j— ^- / / exp (-r 2 + 2rf cos 6)r m drd6 < / (r + f ) m exp (-r 2 )dr 

e m J0 J0 £ m J_ fo 

As f < 1, there exists a constant C(m) > depending only on m such that f[6"9l holds true. 



B.4 Proof of Lemma kL4 



We use polar coordinates and get 

p.v. / — \\3 — d"q = — — p.v. / dr / exp(2rr cos(# — U )) I . . 1 dv 

Jm 2 I V I 5 Jo r Jo \smvj 

e ~ fl r e ~ T \ t a to - a^J cosd ° 
p.v. / dr I cos U exp (2rr cos 6 ) dU 



e 2 Jo r J \sm6 

_f2 

^M 3 Al (f )77o, 
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where we recall that r] = r I cos ^° j anc i where we have set 



sin #o 

oo — r 2 f2TT 



^ poo — r pzir 

M^{z) := - p.v. / dr cos0 exp(2r z cos 9) d6. (182) 

z Jo r Jo 

Standard computations yield that 



i poo ~r z Jf" (n \k k p2n 

4 r 00 e~ r2 , ^ ( 2r )2 P +i^ P +i 

= - P-v- / dr > 7 2p+2 

z -A> r ^ (2p+l)! 



2 2p+1 1 

p=0 



/ 1 7r/2 /'OO 

where ifc := / cos k 9 d6 is the Wallis integral and := / t s ~ 1 e~ t dt is the Gamma 

Jo Jo 

(2p)\ ix 



function. Using the fact that Ii v = , , NO — , we have 

2 2 P(p!)' J 2 



oo 



M 3 Al (z) = 7rV Efe±il ^. (183) 

The radius of convergence of M^ 1 is clearly infinite, since 

r(p+|)(p + l)!(p + 2)! (p+l)(p + 2) 
lim — ^ — 5—— = 1 = oo, 

r(p+|)j?!(p + i)i p + l 

where we have used the standard fact that T(z + 1) = zT(z) for dt(z) > 0. Lemma T4.4I is thus 
established. 



B.5 Proof of Lemma 14.7 

One has 

M*(z) 



poo pzn poo pZn oo lc\ \k k 

\ e~ r2 dr / cos 2 #exp(2r2cos#)c^= / dr / cos 2 gV l / cos fc flrffl 

io Jo Jo Jo k=0 k. 

poo (n \k„k p2tt 00 o2p^2p roo 



fc=0 " u p=0 

00 



00 o2p+l 



Using the fact that I 2p = n0nt n9 T7; we have 



(2p)! 7T 

2 2 ^!) 2 2 
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The radius of convergence of M^ 1 is infinite since 

(2p + 2)(2p+l)((p + 2)!) 2 r(p+i) 
P 4?L (2p + 4)(2p + 3) ((p + l)!) 2 r(p + | + 1) °° ' 



Let M 5 Al (z) be defined by 



/*OC /»27T 

M^(z):= exp(-r 2 )dr / exp(2rz cos 9) d9. (185) 
Jo </o 



We have 



M^(z) = / exp(-r 2 )dr / exp(2rzcos #)d# = / e~ r2 dr V 1 ' \ cos k 9 d9 
Jo Jo Jo k=0 k- Jo 



2 1 2p 

p=0 K t 1 p=0 K ^ 



It is clear that the radius of convergence of M^(-) is infinite. Since M^iz) = M^(z) - 
M^(z), the radius of convergence of M 2 Al (z) is also infinite. 

We now prove that (z) is well-defined and not identically equal to zero. Indeed, 

M^(z)-z 2 Mp{z)-M*\z) 

= 2M*(z)-*±^±p-z**-z*M?\z) 
p=0 W'' 

^ (p + l)(2p + 1) 1 1 2 ^ p 1 2 

3tt ^ P r(p-j) 2p 



2 ^(p + i)(pl)2 z 
Then, the function (z) is defined by 

IV U W- 2 2_ 
which is clearly a non zero entire function. 



2 ^(p+2)((p+i)i; 



186) 



B.6 Proof of Lemma [572] 

We give in this section explicit expressions of ai, 02, and CI3 defined respectively in Eqs. (11281) . 
(11291) . and (I130p . The computations are lengthy but straightforward. 
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We start by computing a\. 
f oc ( 1 

/ -r-^rj T F x r}{tl)(x),T])r}dr} 

roo p2tv 

= ee- f2 °Rg / e- r2 r 2 dr / e 2rf ° cos6 {F£ cos 2 + 2F% cos sin + F 22 (l - cos 2 0)) 
Jo Jo 

I cos 2 6* sin 9 cos #\ T , . 

^sintfcostf l-cos 2 ^^^^ 



/*co /*27t 

ee-^R eo / e _r VcZr / e 2rr "° cos ' 



F^>04 



JO JO 

cos 2 + (F^ 1 - F 22 ) cos 4 (9, 2F e 12 cos 2 9(1 - cos 2 (9) \ 

2F£ cos 2 (9(1 - cos 2 (9), F 22 + (F* 1 - 2F 22 ) cos 2 9 - (ift - F 22 ) cos 4 flj 



The functions M 6 (-), M 7 (-), and M 8 (-) were defined in Eqs. ffT32|) . ffT33]) and flTMft respec- 
tively. Then, we have 

f ^ v T F xV ^(x),v)vdv = ee- f oR eo M(ro)ReMx), (187) 
Jm 2 \V\ 



with 

M(f ) 
._ « 2 



JO) 

< 22 M 7 (f ) + (F^ 1 - F 22 )M 8 (f ), 2F, 12 (M 7 (f ) - M 8 (f )) 

2^ 2 (M 7 (f ) - M 8 (f )), F 22 M 6 (f ) + (F^ 1 - 2F 22 )M 7 (f ) - (F^ 1 - *g»)Af 8 (f ) y 



Then, 

Re M(f )Rl 



Mu+M 2 2 r M u -M 2 2 /cos2fl sin2# \ /-sin2# cos2# N 

2 2+ 2 Vsin2^ - cos 29 ) + 2 \ cos 29 sin2# 0y 



with 



M "+ M » = ^(F 22 M 6 (r- ) + (F--F 22 )M 7 (r- o) ; 

• Mll ~ >l22 = \ (2{Fg - F 22 )M 8 (f ) - (F e \l - 3F 22 )M 7 (f ) - F 22 M 6 (f )) , 
M 12 = 2F e 12 (M 7 (f )-M 8 (f )). 

We also note that 

_ T 1 1 /cos 2# sin 2# \ 

Wo - 2 /2 + 2Vsm^o -cos2fl 0/ J' 

_ ± _ T 1/0 -A l/-sin20 o cos2# N 

- 2 ^ o J + 2 V c °s2^ sin2# 0y 
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We get R 9o M(r )Rj o = M 22 h + (Mu - Mm)wi ~ M i2 + 2M 12 r) L r) r , whi ch 



implies that 

m . i n ( 

-77 F x rj(ip(x),r])ridr] (188) 



On the other hand, one has 

^^-7] T F x r](ifj(x),r])dri 

- e- f2 °ib T (x)R« e~ r2 rdr t e 2rfo cos9 ( F °° °° S 9 + ^ ~ F ^ cosS 6 \ d6 

- e ip{x)K eo J^ e rdr ^ e ^ 2F, 12 cos#(l - cos 2 0) J ^' 

The functions M 9 (-) and M 10 (-) were defined in Eqs. (11351) and (I136P respectively. Then, 
we have 

a e{v) Tn / ,/ n \ 7 =a (T/ N „ fF 22 M 9 (f ) + (F} 1 - F 22 )M 10 (f ) 

Since ip T (x)Rg = ((if}(x),f} ), (^(x),^)), we obtain 

^M^F^^^), 77)^77 = ( [F 22 M 9 (f ) + (F^, 1 - F 22 )M 10 (f )] (^(x), jfo) (189) 

+ 2F e 12 (M 9 (f ) -M 10 (fo))(^(x),^->). 

One also gets 



a£{v \ T F xV d V = ee-'o ^ e~ r \ 2 dr e 2rf0CO8e (F^cos 2 e + F^sm 2 e 

'0 io 



/*oo r2n 

ee- f o / e -Vdr / e 2rf ° cos6 '(F 22 + (F^ 1 - F 22 ) cos 2 
£ e- ? o(F 22 M 6 (f ) + (F£ - F 22 )M 7 (f )). 



Finally, one derives 



^-(rj T F x r])rjdr] 

_ f2 „ f°° _ r 2 . f 2w 2r , cos9 (Ff cos 9 + (F, 11 - F 22 ) cos 3 0\ ^ 
= e r °F 0o / e r rdr e ^r co S e e v e 00 J de 

Jo Jo V 2F^(cos0-cos d 0) y 

Ke °{ 2F e 12 (M 9 (f ) - Mio(fo)) 
Since f^Re^ = (1,0), we have 

/ T#(^)(?/o,^ = ee-'o(F 22 M 9 (f ) + (F^ - F 22 )M 10 (f )). (190) 
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In summary, we get 

at = — — (M 2 2^(x) + {Mn -M 22 )(il>(x),fjo)fjo - Mi^x) 1 + 2Mi 2 (i(f{x) ) fjo)fjo 

[(F 22 M 9 (f ) + (F£ - F, 2 o 2 )M 10 (f ))(^(a;),r /o ) + 2F e 1 () 2 (M 9 (f ) - Mi (f ))^(x),fj^)]fj 

+ (2^(M fl (f ) - M B (f )) + (Fj 1 - F 22 )(M 7 (f ) - M u (f )M*)) 

1 e" f o / 
4ir e \ 

[2F 22 M 6 (f ) + (2^ - 3F 22 )M 7 (f ) - (F^ 1 - F 22 )(M 8 (f ) + M 10 (f )) - F 22 M 9 (f )]^(x) 

- [(F£M 9 (f ) + (f* 1 - F 22 )(M 10 (f ) - 2M 8 (f )) + (F» - 3F 22 )M 7 (f ) 
+F 22 M 6 (f )) (^(x), + 2F e 1 ( f(M 9 (f ) - M w (r ))(i/>(x), f}^)]f} 

- 2F e 12 (M 7 (f ) - M 8 (f ))^(x) ± +4F, 1 o 2 (M 7 (f ) - M 8 (f ))(^(x), ft)^) . 



Let us now compute a 2 . Using the computations performed for the term ai, one has 



|77|3 



rj T F x r]dr] 



-f 2 
e o 



-f 2 
e o 



2tt 



e" r rfr / e 2rf0 cos ^Fi 1 cos 2 + F 22 (l - cos 2 6))d6 



*%M£*{r ) + (F" - F 22 )M^(ro) 



->22\ ■\/ r A- L i- 



The other contribution in a 2 is given by the following expression. 



■ 1 3 ■(ip(x),r])F x 7]dr] 



-f 2 
e o 



-F T R 



e r dr 

M^(ro) 




2- 



3 2rro cos ( 



fl / cos 2 6> sin 9 cos , 

sin 9 cos 8 1 — cos^ / " ( l ' 







i 



R T eMx) 



cos 2# sin 200 
sin 2#o — cos 26* 



-F x (-M^(f )I 2 + (M^(fo) - 2^(r )) ^ 

r 2 

-F x ((M^(r ) - M^(f ))/ 2 + (2M^(f ) - M 5 Al (f ))M T ) 



Therefore, we have 

1 e" f o 



a 2 



47T £ 



+FMM^(f ) - M^(f ))ij(x) + (2M^(f ) - Mf 1 (r )){i/>(x),fjo)fjo 
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Finally, 03 is computed as follows. 

8tt J R 2 \tj\ \v\ 

r>00 /"27T 



1 


-f 2 
e 


8tt 


£ 


1 


-r 2 

e 


8tt 


£ 


1 


-r 2 

e 


8tt 


£ 



ft,, fV*dr / ( . C ° S2,, . .f So "» ) dBRHix) 

J J Vsm0cos0 (1 — cos d) I 

M M^(f )-M^(f )J^ (x) 
(M 5 Al (f ) - M^(f ))^(x) + (2M^(f ) - M^(f ))^(x)^ )% 



This ends the proof of Lemma 15.21 

B.7 Proof of Lemma [573] 

Recall that 

f27T 



Then, one gets 



/ cos 2p 9d9 = AI 2p , / e - r V p rfr = -r(p+-). 
Jo Jo 2 2 



/*oo t>27r 00 r)fc /*oo / i 27r 

M 6 (z) = / e~ r \ 2 dr \ e 2rzcos6 d9 = V -7 [ / e~ r \ k+2 dr cos k 8d6]z 
Jo Jo k=0 k- Jo Jo 



00 o2p+l q 



p=0 





/*OG /*27T 00 o/c /*00 ^277 

M 9 (z) = / e~ r \dr e 2rzcos0 cos 6d6 = V — [ / e~ r \ k+1 dr cos k+l 6dd]z k 

Jo Jo ,_ n k- Jo Jo 



k=0 

00 o2p+2 q 

£(2^ r (p+!> WP+1 



p=0 
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and finally 

"2tt 00 r)k poo p2w 



poo pzn Ok poo pzn 

Afio(z) := / e~ r2 rdr e 2rz cos 9 cos 3 0d0 = V — [ / e~ r \ k+1 dr / cos fc+3 

00 o2p+2 q 

4^ (2p+ 1)! 2 ; 2p+4 



p=0 

Therefore, we obtain 

2M 7 (z) - M 8 (s) - Afio(z) - Mj 41 ^) 



y L—Yip + -)h P+ 2 P P 2 z 2p - y -Y(p + -)/ 2p+4 z 2p+1 . 

^ (2p)! ^t- 2 ;2p+2 2p + 3 ^( 2 p+l)! ^-r 2 ;2p+4 



2M 6 ( 2 ) - 3M 7 (z) + M 8 (z) + M 10 (z) - M 9 (z) - 2M^(z) + 2M^(z) 

°° 92P+1 1 fi n 2 , 1 fi , 7 00 92P+2 o I 

= E ^ + -E(27TT)! r(p + §>'~*T4*~' 
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